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^Airports  for  light  aircraft  must  often  be  built  in  areas  where  base  course 
material  is  inadequate  or  economically  unavailable.  Sand  grids  may  provide  an 
economical  solution  in  these  cases.  Previous  work  with  sand  grids  has  identified 
optimum  grid-cell  geometry  based  on  ultimate  bearing  capacity.  Past  observations 
also  suggest  the  nature  of  sand-grid  behavior  for  small  stresses  and  strains.  From 
these  a  model  was  formulated  for  analyzing  the  performance  of  sand  grids  in 
pavement  systems  and  spot-checked  using  field  data  from  a  full-scale,  accelerated- 
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INTRODUCTION 


BACKGROUND 

In  many  cases,  airports  are  located  on  river  banks  or  in  flood- 
plains  which  do  not  possess  good  subgrade  characteristics.  In  some  in¬ 
stances,  the  subgrade  can  only  be  improved  at  great  expense  by  replacing 
weak  material  with  imported  soil,  or  by  stabilization.  Further,  base 
course  materials  may  also  be  inadequate  or  economically  unavailable.  For 
remote  areas  where  suitable  substitute  soils  cannot  be  obtained  and  where 
an  airport  pavement  is  required,  a  technique  by  which  native  material  can 
be  used  to  support  light  commercial  or  general  aviation  aircraft  needs  to 
be  developed.  The  resulting  criteria  and  methodology  may  provide  greater 
economies  in  the  construction  of  remote  airports  by  utilization  of  cer¬ 
tain  on-site  materials. 

Soil  confinement  systems  have  evolved  from  the  use  of  geotextiles 
in  soil  construction.  The  tensile  strength  of  geotextiles  make  them  a 
logical  complement  to  soils,  which  exhibit  little  or  no  tensile  strength. 
Mass  strength  can  be  increased  by  using  geotextile  cells  to  supply  addi¬ 
tional  confining  stress.  Most  effective  with  granular  materials,  this 
concept  has  been  employed  for  many  years  in  the  form  of  sandbags.  Re¬ 
cently,  the  US  Army  Corps  of  Engineers  has  experimented  with  mechanical 
confineme  it  of  loose  or  weak  cohesionless  materials  to  strengthen  them 
for  use  in  temporary  roads.  Developmental  work  has  been  conducted  using 
grids  forming  interconnected  confining  cells,  open  top  and  bottom  (see 
Figure  1).  Various  materials  have  been  tried,  including  paper,  aluminum, 
and  plastic.  High-density  polyethelene  (HDPE)  has  been  found  to  provide 
the  best  combination  of  strength,  service  life,  and  economy.  The  cell 
system  contributes  to  the  composite,  vertical  compressive  strength  by  in¬ 
creasing  the  horizontal  confining  stress  acting  on  the  cell  filler  mate¬ 
rial.  When  the  filler  material  is  sand,  these  grid  systems  are  sometimes 
known  as  sand-grid  confinement  systems,  or  simply  sand  grids. 

SCOPE 


This  report  presents  a  discussion  of  previous  work  on  confined  soil 
systems  for  pavements  and  the  development  of  a  performance  model.  Previ¬ 
ous  and  concurrent  work  on  soil  confinement  systems,  existing  pavement 
system  performance  models,  and  basic  soil  mechanics  theories  and  prac¬ 
tices  form  the  basis  of  this  analysis.  Laboratory  testing  and  field 
verification  of  the  design  methodology  and  construction  procedure  are  be¬ 
ing  considered  in  a  continuation  of  this  research  effort. 

OBJECTIVE 

The  objective  of  this  report  is  to  develop  an  analysis  of  the  be¬ 
havior  of  soil  confinement  systems,  identify  a  mode  of  behavior  appropri¬ 
ate  to  sand  grids  in  pavement  systems,  and  describe  the  formulation  of  a 
mathematically  based  model  to  predict  sand-grid  pavement  performance 
(serviceability).  Performance  is  to  be  determined  as  a  function  of  air¬ 
craft  landing  gear  configuration,  wheel  loads  and  tire  pressures,  traffic 
intensity,  and  increased  strength  of  the  confined  material. 


CONSTRAINTS 


Pavements  are  traditionally  analyzed  in  layers  according  to  mate¬ 
rial  strength.  Soil  confinement  systems,  because  of  their  moderate 
strength,  may  find  application  in  the  base  course  or  subbase  elements  of 
traditional  pavement  systems,  or  even  as  surface  layers. 

As  noted  earlier,  soil  confinement  systems  are  most  effective  with 
granular  materials.  This  is  because  an  increase  in  strength  with  in¬ 
creased  confining  pressure  is  characteristic  of  granular  materials.  Both 
coarse  and  fine  granular  materials  are  suitable  for  use  in  soil  confine¬ 
ment  systems.  Fine  aggregate  is  often  readily  available  even  when  coarse 
aggregate  is  scarce  or  unavailable.  Hence,  this  study  concentrates  on 
fine  aggregate  or  sand  as  the  cell  filler  material. 

These  material  constraints  are  typical  of  remote  areas  with  rela¬ 
tively  undeveloped  transportation  networks  and  many  light  utility  air¬ 
ports.  Air  traffic  at  these  airports  is  characterized  by  gross  weights 
of  less  than  30,000  lb,*  and  usually  less  than  12,500  lb.  Various  air¬ 
craft  can  be  related  to  existing  pavement  performance  models  and  data 
using  an  equivalent  single-wheel  load  (ESWL)  concept,  such  as  the  one  de- 
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veloped  and  used  by  the  US  Army  Corps  of  Engineers.  The  ESWL  is  ap¬ 
plied  to  the  pavement,  including  the  sand-grid  system,  through  the 
pavement  surfacing  or  wearing  coarse,  if  any.  Light  utility  airport 
pavements  frequently  include  a  2-in.  asphaltic  concrete  wearing  course. 
This  surfacing  serves  to  distribute  both  the  vertical  contact  stress  from 
aircraft  tires  and  the  horizontal  shear  stresses  generated  by  aircraft 
braking  and  turning.  The  base  course  is  thus  protected  from  localized 
material  displacements  or  failures  caused  by  high  tire  contact  pressures 
or  aircraft  braking  and  turning  forces. 

Pavement  performance  is  most  usefully  evaluated  in  terms  of  ser¬ 
viceability.  The  end  of  a  pavement’s  service  life  is  usually  marked  by 
excessive  permanent  deformation  (rutting)  or  cracking.  Brabston,  Barker, 
3 

and  Harvey  found  that  a  strain-based  criteria  provided  the  most  accurate 
and  consistent  prediction  of  flexible  pavement  service  life.  Their  de¬ 
sign  method  is  based  on  limiting  vertical  compressive  strain  at  the  top 
of  the  subgrade  and  horizontal  tensile  strain  at  the  bottom  of  the  bitu¬ 
minous  concrete  layer.  A  review  of  various  strain  criteria  is  given  by 

2 

Barker  and  Brabston. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  ii. 


3 


PREVIOUS  WORK 


DESCRIPTION 

21 

Early  attempts  to  model  the  behavior  of  sand  grids  have  been  ham¬ 
pered  by  a  lack  of  quantitative  data  describing  actual  behavior. 

21 

Mitchell,  Kao,  and  Kavazanjian,  however,  proposed  several  possible 
failure  modes.  These  included  (a)  cell  penetration  into  the  subgrade, 

(b)  cell  rupture,  (c)  cell  wall  buckling,  (d)  simple  bearing  capacity, 

(e)  beam  or  raft  bending,  (f)  deterioration  with  time,  and  (g)  rutting 
under  traffic.  They  then  conducted  a  series  of  model  tests  to  gather  ad¬ 
ditional  quantitative  data.  This  series  formed  an  extension  of  earlier 

26 

work  by  Rea  and  Mitchell.  Like  the  tests  conducted  by  Rea  and 
26 

Mitchell,  these  experiments  were  primarily  concerned  with  the  plate¬ 
bearing  capacity  of  sand  grids  for  various  grid  geometries. 

26  2 1 
Rea  and  Mitchell  and  Mitchell,  Kao,  and  Kavazanjian  suggested 

optimum  geometrical  relationships  based  on  their  experiments.  They  con¬ 
cluded  that  the  width  of  the  grid  cells  should  be  about  two-thirds  the 
width  of  the  loaded  area,  and  the  height  of  the  grid  cells  should  be 
about  four-thirds  the  width  of  the  loaded  area.  These  relationships  are 
consistent  with  classical  soil  mechanics  theory.  According  to  Prandtl's 
1 8 

theory,  the  zone  of  shearing  failure  under  a  loaded  area  extends  to  a 
depth  of  less  than  four-thirds  the  width  of  the  loaded  area  in  soils 
whose  angles  of  internal  friction  are  less  than  about  30  deg.  An  angle 
of  internal  friction  of  30  deg  is  a  typical  design  value  for  sands. 

From  the  standpoint  of  repetitive  loading  due  to  traffic,  the 
plate-bearing  capacity  of  a  sand  grid  can  be  thought  of  as  the  maximum 
load  for  one  repetition.  Thus,  these  optimum  geometrical  relationships 
can  be  used  to  estimate  the  cell  dimensions  required  for  a  particular 
type  of  traffic.  Cell  dimensions  can  then  be  adjusted  based  on  the  per¬ 
formance  of  sand  grids  under  traffic.  After  numerous  tests  using  truck 
30-33 

traffic  Webster*  has  adopted  a  cell  height  of  8  In.  and  a  mean  cell 

diameter  of  7  in.  The  typical  truck  tire  has  a  width  of  about  8  in. 

Thus,  Webster  has  adjusted  his  cell  width  to  approximately  seven-eighths 
the  width  of  the  loaded  area,  and  his  cell  height  to  approximately  equal 
the  width  of  the  loaded  area.  These  dimensions  have  proven  satisfactory 
under  extended  traffic  tests,  and  have  become  the  de  facto  industry  stan¬ 
dard.  They  are  therefore  used  as  a  basis  for  analysis  and  point  of  de¬ 
parture  for  this  study. 

ANALYSIS 

2  6 

Rea  and  Mitchell  noted  that  sand  grid  layers  produce  a  greater 
lateral  distribution  of  vertical  stress  than  unconfined  layers  of  granu¬ 
lar  materials.  This  indicates  that  they  behave  like  beams  or  plates 
which,  unlike  layers  of  unconfined  granular  materials,  can  carry  tensile 


*  S.  L.  Webster.  1983.  Private  Communication,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Miss. 


stresses.  This  observation  suggests  that  the  sand  in  the  unloaded  sand 
grid  is  held  in  compression  by  the  grid  cells  much  like  concrete  is  held 
in  compression  by  the  prestressing  tendons  of  a  prestressed  concrete 
beam.  Then,  under  moderate  loads,  the  composite,  sand-grid  layer  is  able 
to  perform  as  a  flexural  element  without  developing  tensile  strains  in 
the  sand. 

2 1 

Mitchell,  Kao,  and  Kavazanjian  reported  that  the  load-settlement 
behavior  of  plates  on  sand  grids  is  approximately  linear  up  to  moderate 
(50  psi)  stress  levels,  and  that  grid  cells  can  increase  the  effective 
modulus  of  the  sand  filler  material  by  a  factor  of  2  to  3.  They  also 
suggested  a  qualitative  model  for  predicting  the  modulus  of  elasticity  of 
a  sand-grid  layer.  Their  test  results  indicated  that  the  significant  pa¬ 
rameters  are  layer  geometry,  loaded  area-grid  geometry,  sand  modulus, 
grid-material  modulus,  subgrade  modulus,  and  the  number  of  grid  joints 
per  unit  area. 

Hicks ^  concluded  that  the  apparent  Youngs  modulus  of  granular  ma¬ 
terials,  E  as  observed  in  triaxial  tests  can  be  represented  by 
equations 


(1) 


and 


(2) 


where 

kj,  k^,  k^,  k^  =  regression  constants 

o3  =  minor  principal  stress 

1  =  bulk  stress  (sum  of  the  three  principal  stresses) 

Further,  Poisson's  ratio,  v  ,  as  observed  in  triaxial  compression  stress 
states  can  be  described  as 


Parameters  affecting  the  values  of  the  regression  constants  in  these 
equations  include  aggregate  density,  gradation,  soil  type,  degree  of  sat¬ 
uration,  and  particle  angularity.  Thus  the  apparent  E^  and  v  can  be 

estimated  from  the  results  of  laboratory  tests  on  a  particular  material 
prepared  to  simulate  the  anticipated  field  conditions. 

The  problem  is  complicated  by  the  fact  that  stresses  in  pavement 
layers,  during  the  passage  of  traffic  wheel-loads,  vary  widely  in  both 
magnitude  and  orientation.  Stress  conditions  under  passing  wheel-loads 
are  an  area  currently  not  well  understood.  The  reasonably  well  developed 
earth-pressure  theory  applies  to  static  conditions  at  failure. 

The  presence  of  grid  cells  in  the  sand  layer  modifies  the  apparent 
E^  and  v  by  altering  the  stress  state  of  the  soil.  The  combination  of 

the  grid  cells  and  compaction  (by  rollers  during  construction  and  by 
traffic)  apparently  increases  the  horizontal  stress  in  the  sand,  and 
changes  the  shape  of  the  shear  zone.  The  magnitude  of  this  stress  in¬ 
creases  and  the  modified  shape  of  the  shear  zone  is  uncertain.  However, 
the  significant  parameters  influencing  the  process  include  the  elastic 
modulus  of  the  grid  material,  the  geometry  and  dimensions  of  the  grid 
cells  and  the  loaded  area,  the  degree  of  densif ication  of  the  grid 
filler  aggregate,  and  the  intensity  of  the  load.  Preliminary  data  from 
Webster*  (Figure  2)  indicate  that  the  product  of  the  thickness  and  the 
modulus  of  elasticity  of  the  grid  material  is  a  dominant  parameter  influ¬ 
encing  the  rutting  of  a  sand-grid  pavement.  Additional  data  are  required 
to  completely  define  this  relationship,  even  for  a  limited  combination  of 
sand,  traffic,  and  grid  type. 

Finally,  the  apparent,  composite  modulus  of  sand-grid  layers  has 
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been  observed  by  Mitchell,  Kao,  and  Kavazanjian  to  depend  upon  the  mod¬ 
ulus  of  the  material  supporting  the  sand  grid  layer.  This  behavior  can¬ 
not  be  modeled  by  a  methodology  which  considers  only  the  physical  charac¬ 
teristics  of  the  materials  in  the  sand-grid  layer.  Sand-grid  layers  must 
be  analyzed  as  part  of  a  system. 

A  predictive  model  for  sand-grid  moduli  must  consider  three  catego¬ 
ries  of  parameters:  (a)  sand-grid  materials,  (b)  support  conditions,  and 
(c)  loading  conditions.  As  previously  noted,  each  of  these  categories 
consists  of  many  parameters  having  a  significant  impact  on  the  perfor¬ 
mance  of  a  sand-grid  pavement  system. 

EVALUATION 

The  composite  modulus  of  elasticity  of  sand-grid  pavements  can  be 

obtained  from  test  sections  using  nondestructive  testing  (NDT)  methods.^ 
This  has  been  done  using  a  sand-grid  test  section  at  US  Army  Engineer 
Waterways  Experiment  Station  (WES).  The  test  section  is  part  of  a  con¬ 
current  WES  research  project.  It  consists  of  an  8-in. -thick  sand-grid 
layer  over  a  compacted  sand  subgrade. 


*  Webster,  ibid, 
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1  2 
PERMANENT  SURFACE  DEPRESSION,  IN. 


Figure  2.  Effect  of  grid  thickness  and  elastic  modulus  on 
permanent  deformation* 


A  Falling  Weight  Def lectometer  (FWD)  was  used  to  generate  deflec¬ 
tion  data  for  the  BISDEF  layered  elastic  computer  program.  This  propri¬ 
etary  program  is  similar  to  the  CHEVDEF  programs,^  except  that  BISDEF 
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uses  the  BISAR  subroutine  as  the  response  model  to  compute  predicted 
deflections.  Values  for  the  layer  moduli  are  varied  to  match  the  deflec¬ 
tions  predicted  by  the  BISAR  program  to  those  measured  by  the  FWD.  Best 
agreement  between  actual  and  predicted  deflections  was  obtained  using 
modulus  values  of  26,000  psi  for  the  sand-grid  layer  and  16,000  psi  for 
the  subgrade . 


*  Webster,  ibid. 
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APPLICATION 


Webster*  concluded  that  sand  grids  can  provide  soft  subgrades  as 
much  protection  against  the  effects  of  large  vertical  stresses  applied  at 
the  surface  as  unconfined  gravel  layers  up  to  1.6  times  as  thick.  A  sim¬ 
plified  design  procedure  thus  consists  of  substituting  an  8-in.  sand-grid 
layer  for  12  in.  of  gravel  in  any  conventional  design  procedure. 

8 

Coatzee  has  suggested  that  the  8-in.  sand-grid  layer  is  equivalent 
to  at  least  6  in.  of  crushed  aggregate,  based  on  a  layered  elastic  analy¬ 
sis  comparing  sand-grid  pavement  systems  to  conventional  pavement  sys¬ 
tems.  He  checked  the  layered  elastic  stress  predictions  using  a  simpli¬ 
fied  three-dimensional  finite  element  computer  program. 

22 

Nixon  and  Partners  have  adapted  a  design  procedure  based  on  lim¬ 
iting  the  maximum  allowable  vertical  compressive  strain  in  the  granular 
layers  beneath  the  sand  grids.  Stresses,  and  hence  strains,  are  computed 
by  modeling  a  transformed  pavement  section  as  a  semi-infinite,  homoge¬ 
neous,  elastic  half-space.  They  assumed  an  elastic  modulus  of  58,000  psi 
for  the  sand-grid  layer. 

32 

Webster  has  demonstrated  a  construction  procedure  for  temporary 
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access  roads  subject  to  heavy  truck  traffic.  Nixon  and  Partners  have 
outlined  construction  procedures  incorporating  sand  grids  as  base  courses 
or  subbases  under  a  variety  of  pavement  surfacings  and  for  various  types 
of  car  and  truck  traffic.  These  were  used  as  a  basis  for  airport  pave¬ 
ment  construction  procedures  developed  during  this  study. 


*  Webster,  ibid. 
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DESIGN  MODEL  FORMULATION 


PAVEMENT  RESPONSE  MODELS 

g 

Crawford  and  Katona  have  prepared  a  state-of-the-art  report  on  the 
prediction  of  pavement  response.  In  their  discussions,  they  refer  to 
three  types  of  idealizations  of  pavement  structures.  These  are  the 
Westergaard,  layered  elastic,  and  finite  element  idealizations.  To  these 
primary  idealizations  should  be  added  several  significant  mutations  and 
combinations  of  the  three  primary  idealizations. 

Hudson  and  Matlock^  developed  a  model  that  essentially  follows  the 
Westergaard  idealization  but  uses  a  numerical  technique  for  solving  the 
equations  of  bending  for  the  thin  elastic  slab  representing  the  surface 
layer  which  overlies  a  Winkler  foundation.  The  numerical  technique  is 
based  on  finite  difference  approximations  of  continuous  functions,  and 
the  corresponding  physical  idealization  of  the  elastic  pavement  element 
is  similar  to  the  finite  element  idealization.  This  idealization  will  be 
referred  to  as  the  discrete  element  idealization.  A  model  developed  by 

27 

Saxena  combines  the  discrete  element  idealization  of  the  elastic  pave¬ 
ment  element  with  an  elastic  solid  Idealization  of  the  underlying  mate¬ 
rial  rather  than  a  dense  liquid  idealization  (Winkler  foundation).  The 
elastic  solid  idealization  (Boussinesq)  is  a  simplified  version  of  the 

4 

layered  elastic  idealization  in  that  only  one  semi-infinite  layer  is 
considered. 

Huang  and  Wang^  developed  a  model  that  combines  the  finite  element 
idealization  for  thin  elastic  pavement  elements  with  the  dense  liquid 
idealization  for  the  underlying  material.  This  model  has  been  modified 

and  extended  for  Corps  of  Engineer  use  by  Chou^.  A  model  developed  by 

12 

Eberhardt  and  Willmer  is  similar  to  that  developed  by  Huang  and  Wang, 
but  with  an  additional  feature  that  considers  an  intermediate  layer.  A 
procedure  was  developed  in  which  the  top  two  layers  are  modeled  as  an 
equivalent  thin  elastic  plate.  The  finite  element  idealization  is  then 
used  for  the  equivalent  plate  and  the  dense  liquid  idealization  for  the 
remainder  of  the  structure. 

As  stated  previously,  these  four  models  are  simply  mutations  or 
combinations  of  the  three  primary  idealizations  and  are  subject  to  simi¬ 
lar  limitations.  Therefore,  the  following  discussions  are  limited  to  the 
three  primary  idealizations. 

9 

Crawford  and  Katona  provide  detailed  discussions  of  the  three  pri¬ 
mary  idealizations  and  include  discussions  of  various  material  character¬ 
ization  procedures  that  are  necessary  for  quantification  of  properties  of 
the  pavement  structures.  For  the  reader  interested  in  an  in-depth  com- 
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parison,  the  report  by  Crawford  and  Katona  is  recommended.  However,  a 
brief  comparison  follows  which  includes  the  primary  reasons  for  selecting 
a  response  model  based  on  the  layered  elastic  idealization. 


In  1942,  Bunnister  described  a  pavement  system  as  being  composed 
of  an  elastic  layer  overlying  a  semi-infinite  elastic  subgrade.  This 
concept  has  been  adopted  and  refined  by  many  researchers.  For  the  lay¬ 
ered  elastic  idealization  (Figure  3) ,  the  pavement  structure  is  repre¬ 
sented  as  a  series  of  horizontal,  uniform,  elastic  layers  with  properties 

defined  by  (a)  ,  the  modulus  of  elasticity  of  the  iC^  layer; 

(b)  v  ,  the  Poisson's  ratio  of  the  i^  layer,  and  (c)  h  ,  the  thick- 

ness  of  the  i  layer.  These  layers  extend  horizontally  to  infinity,  and 

the  n  layer  also  extends  vertically  to  infinite  depth.  Computer  solu¬ 
tions  can  be  obtained  for  practically  any  multilayer  system.  The  Wester- 
gaard  idealization  (Figure  4)  represents  the  pavement  element  as  a  thin 
elastic  plate  with  properties  defined  by  E^  ,  Vp  *  anc*  ^  ’  over  a 

dense  liquid  (Winkler)  foundation.  The  liquid  foundation  is  character¬ 
ized  as  a  bed  of  springs  having  a  certain  stiffness.  Each  individual 
spring  represents  the  effect  of  the  support  provided  over  a  unit  area. 
This  support  is  quantified  by  a  constant  k  ,  which  is  the  ratio  of  pres¬ 
sure  on  the  unit  area  divided  by  the  deflection.  In  the  basic  Wester- 
gaard  idealization,  loads  were  represented  as  uniform  circular  pressure 
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distributions,  but  procedures  developed  by  Pickett  et  al.  and  Pickett 
25 

and  Ray  permit  pressure  distribution  with  any  shape  to  be  handled. 

The  elastic  layered  idealization  would  appear  to  be  a  more  realis¬ 
tic  representation  of  a  real  pavement  structure  since  pavements  are  truly 
layered  systems,  although  the  materials  may  not  be  truly  elastic.  For 
practical  loadings,  however,  the  materials  can  be  represented  by  quasi- 
elastic  properties.  The  representation  of  the  top  layer  as  a  thin  elas¬ 
tic  plate  (Westergaard  idealization)  or  as  an  elastic  layer  is  equally 
valid  when  the  top  layer  is  a  portland  cement  concrete  (PCC)  slab,  as  in 
rigid  pavements.  The  major  difference  lies  in  the  representation  of  the 
remainder  of  the  structure.  The  use  of  fundamental  constants  E  and  v 
to  represent  the  properties  of  underlying  layers  is  theoretically  more 
sound  than  a  single  constant  k  .  The  Westergaard  idealization  also  sug¬ 
gests  that  there  is  no  lateral  shear  transfer  in  the  underlying  layers, 
an  unrealistic  assumption.  From  a  practical  standpoint,  the  elastic  lay¬ 
ered  idealization  is  also  more  valid.  The  determination  of  k  is  made 
with  a  plate  test  and  represents  the  response  of  the  material  to  a  par¬ 
ticular  loading  condition  (i.e.,  30-in. -diam  plate  and  10-psi  vertical 
pressure),  which  may  be  different  from  that  actually  experienced  in  the 
pavement . 

Experience  has  shown  that,  for  gears  with  closely  spaced  wheels  on 
relatively  thin  slabs  (less  than  about  15  in.),  modeling  of  the  under¬ 
lying  layer  with  a  modulus  of  soil  reaction,  k  ,  produces  reasonable  es¬ 
timates  of  the  response  of  the  pavement.  However,  for  larger  loads 
transmitted  to  the  pavement  through  a  number  of  widely  spaced  wheels,  for 
relatively  thin,  high-strength  (large  stiffness)  base  courses,  and  for 
thick  PCC  slabs,  the  validity  of  the  idealization  decreases.  For  the 
thicker  slabs  and  widely  spaced  wheels,  the  zone  of  influence  (stresses 
in  the  underlying  material)  becomes  much  larger  than  that  under  the 
30-in. -diam  plate  normally  used  in  measurement  of  k  ,  although  the 
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Figure  3.  Elastic  layered  pavement  idealization 


10-psi  contact  pressure  used  may  in  fact  be  valid  for  both  conditions. 
The  effect  of  a  thin,  high-strength  (stiffness)  base  course  will  be  more 
pronounced  on  the  load  deformation  response  of  a  30-in. -diam  plate  than 
on  the  load  deformation  response  of  a  thick  PCC  slab.  The  response  of 
the  30- in. -diam  plate  will  be  significantly  reduced  by  the  thin  base, 
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whereas  the  reduction  in  the  response  of  the  pavement  will  not  be  as 
significant. 

Another  situation  in  which  the  use  of  an  elastic  layered  idealiza¬ 
tion  may  be  more  representative  occurs  when  different  types  of  materials 
at  relatively  shallow  depths  exist  within  the  subgrade  (less  than  20  ft). 
For  instance,  a  stiff  or  a  soft  layer  in  the  subgrade  may  not  signifi¬ 
cantly  affect  the  load  deformation  response  of  a  30-in. -diam  plate,  but 
the  effect  may  be  significant  on  the  load  deformation  response  of  a  thick 
slab  under  a  large  load  on  widely  spaced  wheels. 

Characterization  of  each  layer  with  elastic  constants  obtained  from 
laboratory  tests,  rather  than  one  elastic  constant  obtained  from  field 
tests,  provides  the  designer  greater  flexibility.  Note  that  the  materi¬ 
als  in  pavements  may  behave  neither  elastically  nor  linearly,  but  the 
assumption  of  linear  elasticity  is  made  for  practical  application.  The 
state  of  stress  under  which  the  material  is  tested  in  the  laboratory  may 
be  changed  to  conform  to  the  most  critical  state  of  stress  under  which  it 
may  exist  in  the  pavement.  This  is  contrasted  with  the  constant  state  of 
stress  at  which  the  modulus  of  soil  reaction  is  selected.  There  is  also 
the  flexibility  of  being  able  to  readily  change  the  physical  condition  of 
the  specimen  (moisture,  pore  pressures,  density,  etc.),  whereas  this  can¬ 
not  be  so  easily  accomplished  on  in  situ  material.  Thus,  the  use  of  an 
indirect  correction  for  the  effect  of  saturation  (as  is  done  for  the 

modulus  of  soil  reaction^)  is  not  necessary. 

The  use  of  laboratory  procedures  makes  it  possible  to  test  a  more 
representative  sample  of  the  existing  subgrade  and  a  larger  variety  of 
available  base  course  materials.  With  the  present  design  procedure, 
extensive  plate-bearing  testing  is  uncommon.  Another  factor  to  be  con¬ 
sidered  is  the  repetitive  nature  of  the  loads  applied  to  a  pavement. 
Certainly,  the  use  of  laboratory  tests  will  more  readily  permit  consider¬ 
ation  of  the  effects  of  repeated  load  applications  than  will  use  of  field 
plate-bearing  tests. 

The  assumption  of  completely  bonded  or  completely  frictionless 
layer  interfaces  is  not  considered  to  be  a  significant  weakness.  A  simi¬ 
lar  assumption  is  made  by  the  Westergaard  idealization:  the  interface 
between  the  pavement  element  and  the  underlying  material  is  assumed  to  be 
frictionless.  The  interface  between  a  pavement  element  and  the  second 
layer  is  most  likely  intermediate  between  a  completely  bonded  and  com¬ 
pletely  frictionless  condition.  Between  all  other  layers,  the  assumption 
of  full  bonding  is  probably  more  valid,  being  dependent  on  the  type  mate¬ 
rial  and  construction  procedure.  However,  no  data  exist  to  adequately 
quantify  the  interface  conditions,  although  there  is  a  response  model 
that  can  analytically  consider  intermediate  conditions.  The  computation 
of  the  various  response  parameters  will  certainly  be  affected  by  the  se¬ 
lection  of  the  interface  condition. 

The  comparison  of  the  finite  element  idealization  with  the  Wester¬ 
gaard  and  elastic  layered  idealizations  may  not  be  valid  since  the  finite 
element  idealization  is  basically  a  computation  procedure  rather  than  a 
mathematical  representation  of  the  physical  structure.  As  noted 


previously,  the  finite  element  idealization  may  be  employed  for  the  upper 
layer,  with  a  dense  liquid  or  elastic  solid  representation  for  the  re¬ 
mainder  of  the  structure.  Nevertheless,  consideration  as  a  separate  ide¬ 
alization  has  merit  for  comparisons  between  available  techniques  for  com¬ 
puting  the  response  of  a  pavement  to  loads.  As  discussed  here,  the  fi¬ 
nite  element  representation  considers  the  entire  structure  to  be  broken 
into  a  number  of  finite  elements  (Figure  5). 


NOTE  ELEMENT  NUMBERS  ARE  CIRCLED 
NODE  NUMBERS  ARE  NOT  CIRCLED 
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Figure  5.  Pavement  idealization  using  a  finite  element 
9 

program 

In  the  finite  element  idealization,  the  continuous  pavement  struc¬ 
ture  is  broken  into  a  number  of  elements  connected  at  nodal  points.  The 
material  in  each  element  is  assigned  properties  that  may  vary  from  ele¬ 
ment  to  element.  The  number  of  elements  and  nodal  points  that  may  be 
considered  are  limited  by  computer  capacity,  and  thus  boundary  conditions 
must  be  specified.  The  loads  are  applied  as  concentrated  forces  at  the 
nodal  points.  With  the  aid  of  special  types  of  elements,  discontinui¬ 
ties,  special  interface  conditions,  reinforcing  steel,  and  dowel  bars  may 
be  modeled.  Special  computational  techniques  permit  consideration  of 
voids  and  temperature  and  moisture  gradients  within  the  pavement  struc¬ 
ture.  In  addition,  variable  layer  properties  (thickness  and  load  defor¬ 
mation  properties)  and  nonlinear  material  response  may  also  be  treated. 

However,  there  are  limitations.  For  a  three-dimensional  idealiza¬ 
tion  with  only  a  minimum  number  of  elements  and  refinements  (Figure  5) , 
the  required  time  and  cost  involved  in  applying  the  procedure  to  pavement 
problems  become  prohibitive. 

There  are  plane  strain,  axisymmetric ,  and  prismatic  solid  finite 
element  idealizations,  but  with  all  of  these  idealizations  certain 
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constraints  are  introduced.  If  the  time,  effort,  and  cost  to  apply  the 
models  are  reduced  to  manageable  levels,  the  applicability  to  a  general 
design  procedure  and  improvements  over  simpler  models  are  likewise 
reduced. 

ELASTIC  LAYERED  MODEL 

The  previous  observations  suggest  that  an  elastic  axisymmetric  lay¬ 
ered  model  subjected  to  a  uniform  circular  loaded  surface  area  is  appro¬ 
priate.  This  model  appears  to  describe  the  behavior  of  both  conventional 
pavement  systems  and  subgrade  layers.  Using  a  limiting  strain  criteria 
for  pavement  design  constrains  stresses  and  strains,  determined  by  the 
elastic  layered  analysis,  to  moderate  levels.  Linear  load-settlement  be¬ 
havior  can  thus  be  maintained.  This  model,  which  is  available  in  closed 
form  solutions,  can  also  be  incorporated  into  a  workable  design 
procedure . 

Larger  loads  or  pass  levels  can  be  accommodated  by  adding  a  struc¬ 
tural  surfacing  (e.g.,  more  than  4  in.  of  asphaltic  concrete)  over  the 
sand-grid  layer.  A  three-layer  analysis  is  then  required  to  compute 
strains  and  determine  the  maximum  allowable  aircraft  load.  Such  an  anal¬ 
ysis  is  used  in  the  structural  design  procedure  developed  by  Barker  and 

Brabston. ^ 

Sand-grid  pavements  can  accommodate  much  larger  strains  than  com¬ 
paratively  brittle  bituminous  concrete  pavements  and  still  retain  their 
structural  integrity.  They  are  also  characterized  by  large  (1-  to  2-in.) 
initial  deflections  due  to  densification  of  the  gird  filler  material. 

The  degree  to  which  permanent  deflection  occurs  can  be  controlled  by 
proper  attention  to  compaction.  Since  most  sand  grid  applications  to 
date  have  been  for  military  expedient  low-speed  roads  over  beaches,  lit¬ 
tle  attention  has  been  given  to  achieving  high  degrees  of  compaction  in 
either  the  underlying  layer  or  in  the  sand  within  the  grid  itself.  De¬ 
flections  resulting  from  sand  densification  occur  under  the  first  several 
hundred  traffic  coverages.  Thereafter,  permanent  deflections  accumulate 
much  more  slowly*  (Figure  6).  Thus,  sand-grid  pavements  will  be  struc¬ 
turally  sound  at  levels  of  total,  permanent  surface  deflections  consid¬ 
ered  unsatisfactory  for  bituminous  concrete  pavements.  The  horizontal 
tensile  strain  criteria,  then,  do  not  control  the  performance  of  sand- 
grid  pavements. 

Limiting  the  vertical  compressive  strain  at  the  top  of  the  subgrade 
provides  a  sufficient  failure  criteria  for  sand-grid  pavements.  An  anal- 

2 

ysis  of  the  data  summarized  by  Barker  and  Brabston,  and  shown  in  Fig- 

-3 

ure  7,  suggests  that  limiting  the  resilient  strain  to  0.6  *  10  will  en¬ 
sure  satisfactory  performance  for  up  to  1,000,000  strain  repetitions. 
Alternately,  the  number  of  coverages  to  failure  can  be  included  as  an  ad¬ 
ditional  design  parameter  by  using  one  of  the  relationships  shown  in  Fig¬ 
ure  7,  to  relate  the  calculated  subgrade  strain  to  the  number  of  allow¬ 
able  coverages. 


*  Webster,  ibid. 


TRUCK  PASSES 

(42.000  LB  TAMDEIW  AXLE  TRUCK,  70  PSl  TIRE  PRESSURE) 

Figure  6.  Permanent  deformation  as  a  function  of  traffic 
repetitions* 


The  construction  procedure  must  be  designed  to  prevent  development 
of  unacceptable  surface  roughness  under  traffic.  One  alternative  is  to 
apply  a  large  compactive  effort  during  construction.  Initial  tests  indi¬ 
cate  that  the  sand  filler  material  should  be  compacted  to  about  80  per¬ 
cent  of  its  laboratory  maximum  dry  density.  Another  alternative  is  to 
consider  the  first  several  hundred  traffic  coverages  as  the  last  step  in 
the  construction  and  reducing  surface  roughness,  as  necessary,  by 
applying  a  thin  wearing  surface  process. 


*  Webster,  ibid. 


SAND  GRID  BEHAVIOR 


INTRODUCTION 

As  stated  earlier,  a  model  for  sand  grid  behavior  must  consider 
three  general  features  of  the  problem: 

(a)  Sand-grid  materials 

(b)  Support  conditions  or  subgrade 

(c)  Loading  conditions 

The  model  described  in  this  section  treats  the  sand  grid  and  subgrade 
separately  and  then  compares  the  interface  conditions  for  consistency. 

The  sand  grid  forms  a  composite  structure,  and  the  load  gets  di¬ 
vided  between  the  sand  and  the  plastic  (Figure  8) .  In  the  analysis  of 
the  upper  layer,  both  the  sand  and  plastic  follow  the  deflection  bowl  of 
the  subgrade  surface  computed  from  an  assumed  pressure  distribution.  A 
computer  program  calculates  the  subgrade  pressures  that  act  on  the  sand 
grid  layer  to  deflect  it  in  the  shape  of  the  subgrade  deflection  bowl. 
These  deflections  are  modeled  by  a  piecewise  linear  relationship.  The 
distributed  forces  are  compared  to  the  assumed  pressure  distribution  ap¬ 
plied  to  the  top  of  the  subgrade.  The  calculation  is  then  repeated  until 
the  two  pressure  distributions  are  equal.  The  final  solution  has  both 
equal  vertical  stresses  and  equal  vertical  displacements  acting  at  both 
the  base  of  the  sand  grid  and  the  top  of  the  subgrade  (Figure  8) .  This 
interface  condition  corresponds  to  a  frictionless  or  smooth  interface  in 
a  two-layer  elastic  system. 

The  sand  grid  layer  has  several  characteristics  of  behavior  that 
must  be  addressed  to  successfully  analyze  its  response. 

(a)  The  plan  view  of  the  individual  grid  cells  are  not  simple  plane 
figures  as  shown  in  Figure  8. 

(b)  The  assumed  wheel  load  has  axisymmetry  while  the  sand  grid  ge¬ 
ometry  has  symmetry  about  two  axes. 

(c)  The  sand-plastic  interface  on  the  vertical  sides  of  the  grid 
cells  provide  a  reduced  friction  surface  in  the  layer  making  the  sand 
grid  noncontinuous  above  some  maximum  stress  level. 

(d)  The  repeated  loading  and  unloading  of  the  confined  sand  in  the 
grid  cells  by  compaction  and  traffic  lock  in  lateral  stresses. 

(e)  Since  the  stress  state  in  the  sand  grid  at  some  depth  below  the 
wheel  load  varies  with  the  radial  distance  from  the  application  point, 
the  modulus  of  the  medium  also  varies  in  the  horizontal  plane. 


By  analyzing  the  sand  grid  separately  from  the  subgrade  all  five  of 
these  characteristics  are  treated,  although  with  differing  levels  of 
refinement . 

Figure  8  shows  the  simplified  plan  view  of  the  sand  grid  that  was 
analyzed  with  the  complex  cell  geometry  treated  as  a  series  of  square 
cells  each  measuring  6.4  by  6.4  in.,  standing  8  in.  tall,  and  having  a 
wall  thickness  of  0.055  in.  By  aligning  the  sides  of  the  squares  along 
the  x  and  y  axes  the  behavior  could  be  analyzed  separately  in  the  x  and 
v  directions  and  then  combined  to  describe  the  layer.  Treating  the  8-in. 
deep  sand  grid  layer  as  a  series  of  four  independent  one  dimensional 
( 1 — D )  beams  indicated  as  soil  beam  1,  soil  beam  2,  grid  1,  and  grid  2 
greatly  simplified  calculations  for  this  analysis.  The  assumption  of 
superposition  of  the  solutions  for  orothogonal  beams  is  made  for  ease  of 
solution.  In  reality  the  actual  behavior  is  probably  more  like  that  of  a 
four-way  slab.  Fully  utilizing  symmetry  about  the  x  and  y  axes,  only  the 
x  direction  beams  in  the  positive  x-y  quadrant  were  actually  analyzed. 

The  complete  solution  involved  adding  the  x  direction  beams  to  the  re¬ 
sults  of  the  y  direction  beams  which  were  also  obtained  from  the  x  direc¬ 
tion  calculations. 

This  analysis  did  not  explicitly  include  the  influence  of  the  sand- 
plastic  interface  in  the  calculations.  Applying  a  theoretical  analysis 
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proposed  by  Hadala  to  calculate  the  side  friction  along  the  vertical 
sides  of  a  cylindrical  pressure  chamber  provided  guidelines  for  the  mini¬ 
mum  possible  subgrade  pressure  distribution  directly  below  the  most  heav¬ 
ily  loaded  grid  cell.  This  minimum  pressure  was  less  than  the  value  cal¬ 
culated  from  the  model,  using  values  of  interface  friction  measured  in  a 
direct  shear  test. 

The  stress-strain  behavior  of  sand  is  the  key  feature  in  the  sand- 
grid  response.  The  shear  modulus  of  sand  depends  upon  the  initial  den¬ 
sity,  the  initial  stress  level,  the  strain  level,  and  the  stress  path. 
This  analysis  employs  a  procedure  developed  for  dynamic  analysis  that  de¬ 
termines  a  secant  shear  modulus  during  simple  shear  for  a  particular 
shear  strain  level  and  initial  octahedral  stress  state.  For  the  sand- 
grid  problem,  response  is  treated  as  a  two  step  process;  first  the  wheel 
load  increases  the  initial  stress  level  and  then  the  layer  deforms  under 
the  load.  This  provides  an  approximation  to  the  actual  stress  path  ex¬ 
perienced  by  the  sand,  which  is  different  from  the  stress  path  felt  by 
sand  during  a  simple  shear  loading.  The  octahedral  stress  is  calculated 
using  vertical  stresses  computed  for  the  applied  load  and  the  lateral 
stress  coefficients  computed  from  a  relation  proposed  by  Mayne  and 
20 

Kulhawy,  which  is  described  later.  The  sand  grid  is  initially  analyzed 
using  a  low  strain  shear  modulus.  Following  the  first  calculation  step 
new  moduli  are  computed  from  the  computed  shear  strains  and  compared  to 
original  shear  modulus  values.  The  procedure  iterates  until  the  moduli 
values  used  in  the  program  are  consistent  with  the  computed  shear 
strains . 

Figure  9  shows  a  flow  chart  describing  the  whole  calculation  pro¬ 
cess  involved  in  the  sand-grid  model.  To  start  the  process,  a  pressure 
distribution  acting  on  the  surface  of  the  subgrade  is  assumed. 


ASSUME 

PRESSURE  AS  SIX 
CONCENTRIC  CIRCLES 


STEP  1 


Figure  9.  Flow  chart  for  the  analysis 

Describing  this  distribution  as  a  series  of  six  concentric  constant  pres¬ 
sure  circles,  a  layered  elastic  program  computes  the  subgrade  deflection 
bowl  in  step  2.  (The  results  reported  in  this  report  used  the  computer 
program  RISAR,  since  it  is  considered  a  "benchmark"  for  layered  elastic 
programs.  However,  any  of  several  available  layered  elastic  programs 
could  be  used.)  In  step  3  the  computed  deflection  profile,  the  surface 
load  from  the  applied  wheel  load,  the  resulting  vertical  stresses,  and 
the  values  are  input  into  the  program  SGRID  (Appendix  A) .  The 

weighting  factors  for  load  input  to  the  program  SGRID  are  given  in  Appen¬ 
dix  B.  During  step  4  the  computer  program  SGRID  uses  an  iterative  proce¬ 
dure  to  calculate  the  pressure  distribution  acting  along  the  bottom  of 
the  sand-grid  layer  when  deformed  to  follow  the  subgrade  deflection  ba¬ 
sin.  The  program  is  used  to  analyze  the  response  for  soil  beams  1  and  2 
in  addition  to  ralculating  the  response  of  the  grid  beams  1  and  2  using  a 
simplified  version  of  the  SGRID  program.  The  output  pressures  acting  on 
the  base  of  the  sand-grid  layer  are  compared  in  step  5  to  the  subgrade 
pressure  assumed  in  step  1.  In  step  5  a  simple  computer  program  CAL  (Ap¬ 
pendix  A)  performs  the  process  of  combining  the  x  direction  and  y  direc¬ 
tion  beams  so  that  the  output  pressure  distribution  can  be  compared  to 
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the  input  pressure  on  a  common  basis.  In  step  6  the  two  pressure  distri¬ 
butions  are  compared;  when  input  (step  1)  and  output  (step  5)  stresses 
were  within  ±5  psi,  the  iteration  process  stops.  The  sand-grid  response 
is  then  characterized  by  the  maximum  vertical  subgrade  strain  computed 
during  step  2. 

The  second  general  feature  of  the  problem  is  considered  by  analyz¬ 
ing  the  sand  grid  response  for  three  different  subgrade  modulus  values. 
The  three  solutions  are  called  case  2,  case  3,  and  case  4  representing 
subgrade  modulus  values  of  16,000,  8,000,  and  32,000  psi,  respectively. 
Each  reported  case  is  also  assigned  a  second  number  that  gives  the  itera¬ 
tion  number.  For  example,  case  3-2  refers  to  the  second  iteration  for 
the  situation  with  a  subgrade  modulus  of  8,000  psi. 

The  final  general  feature  of  the  problem,  the  surface  loading  con¬ 
dition,  was  not  varied  for  analyses  reported.  Calculations  assumed  a 
5,600-lb  single-wheel  load  at  100  psi  tire  pressure  acting  on  2-in.  layer 
of  asphalt  overlying  the  sand  grid.  This  load  would  be  typical  for  a 
single-wheel,  light-utility  aircraft  having  a  gross  weight  of  approxi¬ 
mately  12,000  lb.  The  load  on  the  surface  of  the  sand  grids  was  repre¬ 
sented  by  a  10.5  in.  diameter,  65  psi  uniform  load,  and  by  assuming  a 
simple  pyramidal  pressure  distribution  (with  2:1  side  slopes)  in  the  as¬ 
phalt  layer. 

All  of  the  above  calculations  were  performed  on  a  mainframe  com¬ 
puter.  The  computer  capacity  required  to  perform  these  computations  is 
governed  by  the  mathematical  precision  necessary  to  obtain  accurate  solu¬ 
tions  for  the  numerical  methods  employed  by  each  of  the  programs.  The 
minimum  acceptable  capacity  was  not  determined  during  this  study. 

PROPERTIES  OF  THE  COMPONENT  PARTS 

The  plastic  grid  performs  two  functions: 

(a)  It  increases  initial  confining  stress  of  sand  (and  hence,  the 
E  and  G  used  in  the  program  SGRID) . 

(b)  It  carries  some  of  the  load. 

In  this  analysis,  function  1  influenced  the  initial  stress  conditions, 
while  function  2  was  computed  from  the  modified  version  of  the  program 
SGRID.  The  plastic-beam  modulus  values  were  measured  in  the  laboratory. 

The  shear  modulus  was  measured  in  a  torsion  test  according  to  ASTM 
Standard  D  1043.  At  a  temperature  of  70°  F,  the  shear  modulus,  G  , 
equals  32,000  psi.  Using  a  Poisson's  ratio  of  0.5  to  describe  the  plas¬ 
tic  yields  a  stiffness,  E  ,  equal  to  96,000  psi.  When  tested  in  bending 
using  ASTM  D  747,  the  measured  E  value  equals  120,000  psi  for  a  shear 
modulus  of  40,000  psi.  In  analyses  presented  in  this  report,  the  plastic 
is  described  using  the  values  of  G  equal  to  40,000  psi  and  E  values 
equal  to  120,000  psi. 


The  soil  properties  for  the  analyses  presented  in  this  report  rep¬ 
resent  the  concrete  sand  that  filled  sand  grid  cells  in  field  tests 
performed  at  WES.  Figure  10  shows  the  grain  size  distribution. 

Al-Hussaini  and  Perry  *  reported  the  following  properties  for  this  sand: 

Specific  gravity  =  2.66 
Maximum  dry  density  =  117.7  pcf 
Minimum  dry  density  =  98.2  pcf 
Coefficient  of  uniformity  =  2 
Mean  diameter  =  0.5  mm 

The  particles  have  subangular  to  angular  shape.  Nuclear  density  tests 
performed  after  testing  on  a  short  section  of  sand-grid  roadway  measured 
a  dry  density  of  113  to  115  pcf  at  a  water  content  of  4  to  7  percent; 
these  densities  represented  a  relative  density  varying  from  76  to  86  per¬ 
cent.  This  study  used  a  dry  density  of  114  pcf,  for  a  D  equal  to 

81  percent.  To  analyze  the  sand-grid  beam,  the  shear  modulus,  Poisson's 
ratio  and  sand-plastic  friction  angle  must  be  determined.  The  shear  mod¬ 
ulus  and  Poisson's  ratio  were  based  upon  test  results  and  proposed  rela- 
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tions  reported  in  the  literature  '  while  the  sand-plastic  friction  was 
measured  in  a  direct  shear  test. 

The  shear  modulus  for  a  sand  depends  upon  both  the  initial  stress 
state  and  upon  the  shear  strain  magnitude.  These  Influences  can  be 
treated  separately  using  a  stress  dependent  maximum  shear  modulus  and 

14 

then  decreasing  it  consistent  with  the  shear  strain.  Hardin  and  Black 
found  laboratory  shear  moduli  measured  at  low  shear  strains  were  fit  by 
the  following  equation: 


G 

max 


1230 


(2.973  -  e)' 
1  +  e 


(°oct) 


1/2 


(4) 


where 


G  =  shear  modulus,  psi 
max 

e  =  void  ratio,  psi 


o  =  octahedral  normal  stress,  psi 
oct 

For  the  concrete  sand  used  in  this  study,  which  has  a  void  ratio  equal  to 
0.45,  this  expression  becomes: 
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Figure  10.  Grain-size  distribution 
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Seed  and  Idriss  used  equation  6  to  give  G  and  °oct  *n  psf: 


G  =  1000  K_(o  J 
2  oct 


The  constant  K„  is  a  function  of  both  shear  strain  and  relative  den- 
*■ 

sity.  At  the  small  shear  strain  value  of  10  in. /in.,  equals  64 

for  a  sand  having  a  relative  density  equals  to  80  percent.  If  both  G 

and  o  are  given  in  psi,  equation  6  becomes: 
oct 


G  -  5300  (o  V 
max  oct 


Preliminary  computations  were  made  using  equation  8: 

1/2 

G  -  5500  (o  J)  ' 
max  oct 

to  model  the  higher  relative  density  expected  under  aircraft  loads. 


Is : 


The  initial  stress  required  to  compute  G 

max 


oct 
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where  K  is  the  ratio  of  horizontal  to  vertical  stress.  The  vertical 
o 

stress  at  a  point  equals  the  initial  vertical  stress  caused  by  the  over- 
lying  soil  and  asphalt  plus  the  added  stress  from  the  applied  wheel  load. 
20 

Mayne  and  Kulhawy  empirically  related  the  value  of  Kq  in 

soils  to  the  friction  angle,  <f>  ,  and  the  stress  history  by  equation  10. 


K  =  (1  -  sin  0) 
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OCR  (1  - 
max 
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The  quantities  OCR  and  OCR  refer  to  the  stress  history.  The  OCR 

max 

value  equals  the  ratio  of  the  maximum  vertical  stress  to  the  current  ver¬ 
tical  stress.  For  the  sand  filling  the  sand  grid  cells  the  maximum  ver¬ 
tical  stress  is  assumed  to  have  occurred  during  construction,  when  a 

5,600  lb  wheel  load  compacted  the  soil.  The  OCR  equals  the  ratio  of 

max 

the  maximum  vertical  stress  to  the  minimum  vertical  stress.  During  the 
first  unloading  cycle  for  a  soil,  OCR  equals  0CRmax  while  for  later 

load/unload  cycles  the  OCR  value  remains  constant  at  about  40  while 

max 

the  OCR  value  varies.  Figure  11  compares  the  results  of  equation  10  to 
the  measured  horizontal  stresses  during  a  load-unload-reload  test  on  fil- 

20 

ter  sand  reported  by  Mayne  and  Kulhawy. 

The  response  of  soil  in  simple  shear  is  nonlinear  when  shear 

strains  exceed  10  ^  in. /in.,  with  the  tangent  shear  modulus  decreasing. 
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Seed  and  Idriss  developed  the  relationship  shown  in  Figure  12  to  de¬ 
scribe  the  secant  shear  modulus  as  a  function  of  shear  strain.  For  these 
analyses,  the  average  curve,  shown  as  a  dashed  line  in  Figure  12,  was 

used  to  find  the  ratio  G/G  as  a  function  of  shear  strain,  y  .  Corn- 

max 

bining  G/G^y  with  equation  8,  the  shear  modulus  for  the  sand  is  given 
in  equation  11. 
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For  calculations,  the  Poisson's  ratio,  v  ,  equals  0.35.  Table  1  summa¬ 
rizes  the  material  input  properties  used  in  each  of  the  three  analyzed 
cases  or  references  the  appropriate  equation. 
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Figure  11. 


The  analysis  for  sand  grid  response  uses  an  elastic  approach  and 
assumes  no  developed  plastic  failure  mechanisms.  This  holds  true  when 
the  maximum  vertical  shear  stress  does  not  exceed  the  frictional  strength 
between  the  sand  and  the  plastic  grid.  Therefore,  the  friction  of  the 
interface  was  measured  in  direct  shear  tests.  Figure  13  shows  the 
strength  envelope  for  sliding  of  sand  on  sections  of  plastic  grid  trimmed 
from  a  completed  field  test,  and  compares  it  to  the  strength  envelope  for 
sand.  For  plastic-sand  sliding,  the  angle  of  friction,  6  ,  equals 
21  deg,  while  for  sand,  the  angle  of  internal  friction,  $  ,  equals 
40  deg.  The  ratio  tan  6/tan  if)  equals  0.46.  The  plastic  significantly 
reduces  the  sand  layer's  capacity  to  transmit  shear  stresses  along  verti¬ 
cal  planes.  In  the  next  section  the  significance  of  this  reduced  fric¬ 
tion  surface  upon  subgrade  pressure  distributions  will  be  examined. 


SUBGRADE  ANALYSIS 


The  deflection  profile  of  the  subgrade  was  computed  for  an  assumed 
pressure  distribution.  The  initial  pressure  distributions  used  in  the 
analyses  were  based  upon  results  computed  for  a  two-layer  system  with  an 
8  in.  layer  overlying  a  subgrade.  The  area  under  the  assumed  pressure 
distribution  equals  the  total  wheel  load  of  5,600  lb. 


The  pressure  distribution  underlying  the  sand  grid  cell  loaded  by 
the  wheel  is  controlled  by  the  frictional  force  acting  along  the  vertical 
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sand  plastic  interface.  Hadala  presented  an  analysis  to  compute  the 
vertical  stress  acting  at  the  base  of  a  sand  sample  confined  by  a  cylin¬ 
drical  chamber.  Modifying  his  analysis  for  the  rectangular  prismatic 
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Figure  13.  Direct  shear  test  results 


confining  chamber  formed  by  the  simplified  version  of  the  sand  grid  cells 
leads  to  equation  12. 

a  =o  exp  (-4K  tan  6  H/B)  (12) 

v  o 

where 

=  applied  surface  load 

K  =  ratio  between  the  horizontal  and  vertical  stress 

tan  6  =  the  coefficient  of  friction  between  the  sand  and  plastic 
grid 

Figure  14  shows  the  value  of  at  the  base  of  the  sand  grid  cell 

loaded  by  a  65  psi  wheel  load  plotted  as  a  function  of  K  .  For  the  min¬ 
imum  value  of  K  used  in  the  analysis  of  0.38,  the  minimum  subgrade 
pressure  equals  31.3  psi.  For  the  three  cases  evaluated,  the  central 
subgrade  pressure  always  exceeds  34  psi.  Therefore  the  subgrade  values 
exceeded  the  minimum  allowable  values.  Since  other  grid  cells  carried 
small  surface  pressures,  the  minimum  limit  only  applies  for  the  pressures 
underlying  the  central  cell. 

The  subgrade  deflection  profiles  were  computed  using  the  computer 
program  BISAR  although  any  layered  elastic  program  could  be  used.  As 
noted  earlier,  BISAR  was  chosen  for  this  study  because  it  is  considered  a 
"benchmark"  for  layered  elastic  programs.  Figure  15  shows  a  profile  of 
the  geometry  analyzed.  A  rigid  base  under  the  subgrade  was  located  at  a 
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20  ft  depth  as  recommended  by  Barker  and  Brabston.  The  subgrade  pres¬ 
sure  distribution  was  approximated  by  six  concentric  circular  pressure 
rings  to  satisfy  the  input  limitations  of  the  program  BISAR,  and  this 
procedure  led  to  the  step  shaped  pressure  pattern  pictured  in  Figure  15. 

Figure  16  shows  the  computed  deflection  profile  used  for  one  of  the 
analyzed  cases.  Output  values  correspond  to  the  required  input  locations 
for  the  program  SGRID,  described  in  the  next  section.  Displacements  were 
computed  at  points  along  all  four  of  the  beams  that  were  analyzed.  The 
results  of  the  subgrade  analysis  reflect  the  fact  that  the  pressure  was 
step  shaped  instead  of  continuous.  The  calculated  shear  strains  reach 
their  maximum  values  at  x  values  representing  the  edges  of  load  steps. 

SAND  GRID  ANALYSIS 

For  analysis  the  sand-grid  layer  was  divided  into  a  series  of 
1-D  beams.  Each  of  the  beams  was  in  turn  subdivided  into  a  series  of 
subcells.  The  displacements  and  forces  were  calculated  at  only  the  cor¬ 
ner  points  of  the  subcells.  Figure  17  shows  a  schematic  of  the  analyzed 
geometry.  Fully  utilizing  the  symmetry  of  the  problem,  only  the  positive 
quadrant  of  the  x-y  plane  was  analyzed.  The  responses  of  the  four  dif¬ 
ferent  beams  in  the  x  direction  were  calculated  and  doubled  to  get  the 
full  solution  for  the  contributions  from  both  the  x  and  y  directions. 
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PROFILE 


Figure  15.  Applied  pressure  distribution 
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Figure  16.  Deflection  profile 


cross-hatched  area  shown  in  Figure  17  because  they  contribute  an  insig¬ 
nificant  proportion  of  the  total  load. 

The  cross  section  shown  in  Figure  17  portrays  the  geometry  of  the 
analyzed  beams.  Each  grid  cell  was  subdivided  into  four  subcells  with 
the  response  computed  at  the  subcell  corner  points.  The  soil  beams  were 
analyzed  assuming  plane  strain  behavior  while  plane  stress  applied  for 
the  grid  beams.  Distributed  loads  were  converted  to  point  loads  acting 
at  the  subcell  corners.  Each  subcell  corner  carried  the  distributed  load 
applied  over  a  region  extending  from  the  midpoint  between  that  corner  and 
the  preceding  corner  to  the  midpoint  between  that  corner  and  the  follow¬ 
ing  corner. 

Figure  18  shows  the  flow  chart  i.or  the  computer  program  SGRID  used 
to  model  the  response  of  each  soil  beam.  In  step  1,  two  sets  of  informa¬ 
tion  are  input:  first,  the  data  needed  to  characterize  the  response  of 
the  soil,  and  second,  the  prescribed  corner  forces  imposed  by  the  wheel 
load  and  the  corner  deflections  computed  from  the  subgrade  analysis. 
Figure  19  shows  the  input  and  output  data  for  the  SGRID  computer  program. 
Figure  19a  shows  the  information  input  for  the  analysis  of  soil  beam  1  in 

case  4.  For  each  sand  subcell,  the  K  value,  the  vertical  stress 

o 

level,  and  the  Poisson's  ratio  are  required.  In  step  2  the  low-strain 
shear  modulus  values  are  computed  for  each  subcell  using  the  input  values 
in  equation  8.  In  step  3  the  stiffness  matrix  [K]  is  computed  from  the 
geometry  of  the  problem  and  from  the  material  properties.  Using  the 
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methods  described  by  Smith  the  subcell  corner  forces  { F }  can  be  re¬ 
lated  to  the  subcell  corner  displacements  {U}  in  equation  13. 

(F>  =  [K]  {U}  (13) 

In  step  4  the  appropriate,  known  displacement  and  force  values  are  in¬ 
cluded  in  equation  13,  and  the  remaining  unknown  displacement  and  force 
values  are  calculated.  The  subcell  shear  strains  are  then  computed  from 
the  calculated  values  of  displacements  and  forces.  Using  these  shear 
strain  values  in  step  5,  new  values  for  shear  modulus  values,  G^+^  »  are 

computed  from  equation  11.  In  step  6,  the  new  G  values  are  compared  to 
the  values,  ,  used  to  construct  the  [K]  matrix  in  step  3.  If  the 

G^+j  values  are  more  than  10  percent  different  from  the  G^  values,  the 

program  performs  steps  3  through  6  using  the  G^+j  values.  Iteration 

proceeds  until  the  two  G  values  are  less  than  10  percent  different,  at 
which  time  the  results  are  output.  Figure  19b  shows  the  final  results 
for  beam  1  in  case  4-3.  The  values  for  all  three  cases  are  portrayed  in 
Appendix  C. 

ITERATION  AND  CONVERGENCE 

The  SGRID  output  loads  must  be  converted  to  a  pressure  distribution 
before  comparing  them  to  the  input  pressure  distribution.  The  load  for 
each  subcell  corner  along  the  base  of  the  soil  beams  acts  over  an  area 


Figure  18.  SGRID  flow  chart 
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whose  width  equals  the  beam  width  and  whose  length  equals  the  distance 
from  the  preceding  midpoint  to  the  following  midpoint.  The  load  calcu¬ 
lated  for  the  plastic  grid  beam  subcell  corners  was  distributed  only 
along  the  length  from  the  preceding  midpoint  to  the  following  midpoint, 
assuming  negligible  width.  To  compare  these  rectangular  and  linear  re¬ 
gions  of  constant  stress  to  the  ring  shaped  regions  of  constant  stress 
input  to  BISAR,  the  SGRID  results  were  transformed  to  square  regions  se¬ 
lected  to  correspond  with  the  rings.  The  computer  program  CAL  combined 
the  subcell  corner  loads  using  weighting  factors  calculated  from  the  por¬ 
tion  of  the  acting  subcell  rectangular  area  lying  within  each  square  re¬ 
gion.  The  square  regions  making  up  each  ring  were  then  averaged  to  ob¬ 
tain  the  output  pressure  distribution  for  that  ring.  Tables  2,  3,  and  4 
summarize  the  input  and  output  pressures  for  each  iteration  of  cases  2, 

3 ,  and  4 . 

The  iteration  procedure  involved  selecting  the  average  of  the  input 
and  the  output  pressures  to  be  used  for  the  input  pressure  to  the  next 
iteration.  Before  these  values  could  be  used,  the  total  area  under  the 
pressure  distribution  was  adjusted  to  give  an  area  equal  to  the  5,600  lb 
wheel  load.  Figure  20  graphically  portrays  the  input  and  output  pres¬ 
sures  used  in  case  4,  iteration  1.  The  output  values  for  the  center  two 
rings  are  high  by  15  and  10  psi,  respectively.  Horizontal  arrows  show 
the  pressure  values  selected  for  the  second  iteration.  Figure  21  com¬ 
pares  the  input  pressures  to  the  output  pressures  for  iteration  2.  Fi¬ 
nally,  Figure  22  compares  the  pressures  for  iteration  3.  Note  that  the 
output  values  fall  above  the  input  values  by  less  than  4  psi. 

The  key  decision  during  an  iteration  process  involves  the  level  of 
agreement  between  the  input  and  output  values.  For  this  analysis,  itera¬ 
tion  continues  until  the  input  and  output  pressures  agree  to  within 
±5  psi.  Early  experience  showed  both  that  it  is  difficult  to  get  better 
agreement  criteria  than  this,  given  the  scale  of  discretization  for  the 
model,  and  that  additional  refinement  does  not  significantly  influence 
the  maximum  vertical  subgrade  strain  computed  from  BISAR, 

RESULTS 

Table  5  summarizes  the  results  of  the  final  iteration  for  each  of 
the  cases  analyzed.  The  maximum  vertical  subgrade  strain  represents  the 
feature  of  performance  recommended  by  Barker  and  Brabston  for  design  of 
flexible  airport  pavements.  Their  design  plot  showing  the  repetitions  of 
strain,  as  a  function  of  maximum  vertical  strain  at  the  top  of  the  sub¬ 
grade,  was  presented  in  an  earlier  section.  The  design  discussion  in 
this  section  uses  the  lower  bound  to  the  range  proposed  by  Brabston, 

3 

Barker,  and  Harvey  and  extrapolates  it  to  low  repetitions  of  strain. 
Table  5  shows  that  the  allowable  number  of  passes  by  a  5,600  lb  wheel 
(repetitions  of  strain)  varies  from  70  to  400,000  as  the  subgrade  modulus 
increases  from  8,000  to  32,000  psi. 

The  calculated  values  of  subgrade  strain  versus  the  subgrade  modu¬ 
lus  value  have  been  plotted  in  Figure  23.  This  plot  p  ovides  the  asis 
of  a  design  plot  for  a  single  wheel  lo  d  acting  on  an  8  in.  thi  k  sand- 
grid  layer  having  an  average  cell  width  equal  to  6.4  in.  The  vertical 
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Results  for  Case 


CASE  4-1 


RADIAL  DISTANCE.  X.  INCHES 


Table  5 


Summary  of  Results 


Case 

Subgrade  Modulus 
psi 

Maximum  Subgrade  Strain 
in. /in. 

Number  of  Coverages 

2 

16,000 

0.95  *  10-3 

7,750 

3 

8 , 000 

1.7  x  10-3 

71 

4 

32,000 

0.58  x  io"3 

415,000 

subgrade  strain  axis  can  also  be  shown  as  the  number  of  passes  by  using 

3 

the  Brabston,  Barker,  and  Harvey  relationship  to  relate  subgrade  strain 
to  low  repetitions.  This  is  an  appropriate  relationship,  since  it  is 
also  based  on  an  elastic  model  for  subgrade  behavior.  Additional  dashed 
lines  have  been  added  to  Figure  23  to  represent,  conceptually,  the  single 
wheel  load  as  a  third  design  parameter.  These  lines  were  positioned  ap¬ 
proximately,  by  scaling  the  strains  computed  for  a  5,600  lb  wheel  load, 
by  the  ratio  of  the  wheel  loads.  Many  additional  analyses  of  the  type 
described  in  this  report  are  required  to  complete  the  design  curves. 
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Figure  23.  Conceptual  design  chart 


FIELD  VALIDATION 


FIELD  DATA 

Two  full-scale  test  sections  were  used  to  generate  data  for  the  an¬ 
alytical  model  described  in  the  previous  section.  The  first  was  an  un¬ 
surfaced  test  section  constructed  under  a  previous  WES  project.  This 
test  section  provided  material  for  laboratory  and  field  testing  of  the 
physical  properties  of  the  sand,  grid,  and  subgrade.  These  properties 
are  those  described  and  summarized  in  the  previous  section.  The  second 
test  section  was  part  of  a  larger  test  section  which  was  built  as  a  joint 
project  to  provide  field  validation  of  the  analytical  model  under  this 
project  and  to  study  the  performance  of  various  surfacings  and  grid-fill 
materials  under  a  separate  WES  project. 

The  general  layout  of  the  second  test  section  is  shown  in  Fig¬ 
ure  24.  Item  3  was  constructed  to  match  the  input  parameters  used  in  the 
hypothetical  case  studies  described  in  the  previous  section.  Items  4  and 
12  are  similar.  Item  4  has  a  thicker  asphalt  wearing  course,  which 
changes  the  ESWL  on  the  top  of  the  sand  grid.  Item  12  has  a  different 
type  of  sand  for  both  grid  fill  and  subgrade,  which  changes  both  the  sub¬ 
grade  modulus  and  the  distribution  of  stresses  (k-values)  in  the  sand- 
grid  layer.  Although  they  were  built  primarily  to  address  the  objectives 
of  the  other  WES  project,  they  could  provide  data  for  additional,  future 
case  studies  and  refinement  of  the  analytical  model.  Considerable  expen¬ 
ditures  of  additional  engineering  and  computer  time  would  be  required  to 
develop  the  conceptual  design  charts  for  comparison  with  these  test 
items . 


Item  3  consisted  of  2  in.  of  asphalt  wearing  course  over  an  8-in. 
sand-grid  layer  placed  on  a  compacted  sand  subgrade.  Asphalt  thickness 
was  checked  by  taking  core  samples.  Actual  thicknesses  ranged  from  1.5 
to  2.25  in.  The  dimensions  and  properties  of  the  HDPE  grids  are  as  shown 
in  Figure  1  and  described  in  the  previous  section,  except  that  the  thick¬ 
ness  of  the  plastic  was  0.05  in.  for  the  test  section.  As  production 
tolerances  of  about  ±0.005  in.  are  considered  reasonable,  this  difference 
will  be  ignored  for  the  purposes  of  the  analysis  to  follow.  The  sand 
used  to  build  the  subgrade  and  to  fill  the  sand-grid  cells  is  from  the 
same  source  as  that  described  in  the  previous  section.  Its  grain-size 
distribution  is  shown  in  Figure  10.  The  sand  was  compacted  in  the  grid 
cells  to  a  dry  density  of  105  pcf.  The  moisture  content  ranged  from  4  to 
7  percent.  The  subgrade  stiffness  was  evaluated  using  a  FWD  after  con¬ 
struction  was  completed.  The  initial  subgrade  modulus  was  calculated  to 
be  14,000  psi,  densifying  and  stiffening  to  the  expected  16,000  psi  after 
the  first  500  passes  of  the  test  vehicle. 

The  test  vehicle  was  a  standard  5-ton,  tandem-axle,  M51  military 
dump  truck  carrying  the  maximum  highway  payload  of  20,000  lb.  The  re¬ 
sulting  gross  vehicle  weight  was  42,000  lb.  The  truck  was  fitted  with 
standard,  military,  off-road  tires  inflated  to  the  maximum  cold  inflation 
pressure  of  70  psi. 
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SECTION  PERTAINING 


Lgure  24.  General  layout,  sand-erid  test  section 


Natural  vehicle  wander  was  allowed,  and  a  reversal  of  vehicle  di¬ 
rection  every  100  passes  was  used  to  ensure  a  more  symmetrical  pattern  of 
vehicle  wander. 

Cross-section  data  were  collected  at  0,  100,  500,  1,000,  2,000, 
3,000,  4,000,  and  5,000  passes,  and  supplemented  by  photographs  and  FWD 
measurements.  Damage  to  pavement  systems  is  generally  a  function  of  the 
logarithm  (base  10)  of  traffic  intensity.  A  great  deal  of  additional 
trafficking  would  therefore  be  required  to  produce  significant  additional 
data  beyond  that  obtained  at  5,000  passes.  Hence,  traffic  was  halted  at 
5,000  passes. 

No  destructive  field  testing  was  done,  and  additional  trafficking 
by  other  vehicles  is  anticipated  under  the  other  project.  The  average, 
permanent,  differential  surface  deflection  or  rut  depth  was  obtained  from 
the  cross-section  data.  It  is  plotted  as  a  function  of  vehicle  passes  in 
Figure  25. 


Figure  25.  Average  rut  depth  as  a  function  of  vehicle  passes 


ANALYSIS 

The  data  presented  in  Figure  25  suggest  that  initial  consolidation 
or  shakedown  has  occurred  with  the  rut  depth  stabilizing  at  0.3  in.  A 
pavement  failure  (rut  depth  £1.0  in.)  has  not  occurred  and  does  not  ap¬ 
pear  imminent.  This  information  can  be  used  to  check  the  conceptual  de¬ 
sign  chart  shown  in  Figure  23,  if  the  test-section  traffic  can  be  ex¬ 
pressed  in  terms  of  coverages  of  a  single  wheel  load. 

The  effect  of  the  test  vehicle  on  the  test  section  can  be  described 
in  terms  of  ESWL  as  noted  earlier.  Based  on  the  concepts  presented  by 

an  ESWL  curve  for  the  test  vehicle's  18,000  lb  axle,  as  shown 


Pereira , 


in  Figure  26,  can  be  constructed.  From  Figure  26,  the  ESWL  is  seen  to 
vary  from  31  percent  of  5,600  lb  at  the  top  of  the  sand-grid  layer  to 
42  percent  or  7,600  lb  at  the  bottom  of  the  sand-grid  layer.  Thus,  the 
effective  load  induced  by  one  18,000  lb  axle  of  the  test  vehicle  is 
equivalent  to  the  5,600  lb  single-wheel  load  assumed  for  the  development 
of  the  analytical  model,  at  the  top  of  the  sand-grid  layer.  The  effec¬ 
tive  load  is  somewhat  larger  at  greater  depths.  The  test  vehicle  thus 
provides  a  slightly  more  severe  validation  load  than  a  single  5,600  lb 
wheel.  It  can  be  used  to  provide  a  conservative  validation  check  of  the 
analytical  model  developed  for  a  5,600  lb  single  wheel  load. 

Because  the  test  vehicle  has  two  18,000  lb  axles,  and  because  the 
test  vehicle  exhibits  a  lateral  wander  during  trafficking,  one  vehicle 
pass  does  not  provide  at  least  one  coverage,  or  strain  repetition  across 
the  entire  width  of  the  test  item.  A  test  vehicle  "pass-to-coverage"  ra¬ 
tio  is  required  to  compare  the  vehicle  passes  shown  in  Figure  24  to  the 
coverages  used  in  Figures  7  and  23.  Using  vehicle  wander  concepts  pre- 
23 

sented  by  Pereira  a  pass-to-coverage  ratio  for  the  18,000  lb  axle  has 
been  calculated  to  be  2.64.  Since  the  test  vehicle  has  two  18,000  lb  ax¬ 
les,  and  the  effect  of  the  front  axle  can  be  ignored  because  of  its  rela¬ 
tively  small  load,  the  test  vehicle  pass-to-coverage  ratio  can  be  taken 
as  2.64  t  2  =  1.32  . 

Based  on  the  above.  Item  3  of  the  test  section  was  subjected  to 
more  than  5,000  passes  *  1.32  per  coverage  to  yield  3,788  coverages  of 
a  5,600  lb  equivalent  single  wheel  load  without  experiencing  failure. 
Figure  27  shows  this  point  superimposed  on  a  portion  of  Figure  23.  An 
arrow  through  the  data  point  indicates  the  direction  the  point  would 
translate  had  trafficking  been  carried  to  failure.  For  a  subgrade  modu¬ 
lus  of  16,000  psi,  the  analytical  model  predicts  satisfactory  performance 
to  7,750  coverages.  These  two  values  are  quite  close  in  terms  of  the 
logarithm  (base  10)  of  coverages  commonly  used  in  pavement  design.  The 
analytical  model  predicted  satisfactory  performance  to  log  (7,750)  - 
3.9  ,  and  the  test  section  demonstrated  satisfactory  performance  to  at 
least  log  (3,788)  =  3.6  .  It  can  be  inferred  from  these  results  that 
the  model  is  conservative. 

These  promising  results  support  spending  the  time  and  money  neces¬ 
sary  to  construct  additional  test  sections,  having  various  subgrade 
strengths,  and  trafficking  them  to  failure  with  a  variety  of  equivalent 
single  wheel  loads.  Should  additional  testing  indicate  that  the  model  is 
too  conservative,  the  model  is  easily  calibrated  by  selecting  a  less 
conservative  relationship  between  vertical  subgrade  strain  and  allowable 
coverages. 
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and  predicted  performance 


CONCLUSION  AND  RECOMMENDATIONS 


The  model  described  in  this  report  provides  a  means  of  considering  five 
important  characteristics  of  sand  grids.  These  are:  the  complex  geometry,  the 
different  symmetry  for  the  load  than  for  the  grid,  the  vertical  sand-plastic 
interface,  the  importance  of  locked-in  lateral  3and  stresses,  and  the 
variation  of  sand  modulus  in  the  horizontal  plane.  The  model  handles  these 
characteristics  with  differing  levels  of  sophistication.  Although  the  results 
reflect  the  required  simplifications,  they  provide  a  design  result  which 
compares  favorably  with  measured  test  section  behavior. 

The  flow  chart  described  in  the  text,  and  shown  in  Figure  9,  was  used  to 
generate  a  conceptual  design  chart.  Physical  properties  and  dimensions  for 
the  sand  grid  were  taken  from  the  commercially  available  product  that  is  being 
used  to  construct  sand  grid  test  sections  at  WES.  The  material  properties  for 
the  sand  match  those  of  the  concrete  sand  described  in  the  text  and  Figure 
10.  Loads  on  the  sand  grid  were  based  on  a  56C0  lb.  single  wheel  load  acting 
through  a  2  in.  asphalt  concrete  wearing  surface.  This  conceptual  design 
chart  is  shown  in  Figure  23,  including  estimated  lines  for  3000  and  9000  lb. 
single  wheel  loads.  The  number  of  allowable  coverages  shown  on  the  right-hand 


block  of  the  flow  chart. 

A  design  chart  similar  to  Figure  23  can  be  generated  for  any  particular 
combination  of  sand  grid  and  material  used  to  fill  the  grid  cells.  Sand  grid 
pavement  systems  can  be  designed  for  multiple-wheel  gear  configurations  by 
using  equivalent  single  wheel  load  theory  to  determine  an  equivalent  single 
wheel  load  to  use  with  the  single  wheel  load  design  chart.  The  appropriate 
chart  is  then  entered  with  the  known  design  parameters,  to  determine  the  value 
of  the  unknown  parameter  or  acceptable  combinations  of  the  design 
parameters. 

For  example,  suppose  a  light  utility  airport  pavement  is  to  be 
constructed  using  the  materials  used  to  develop  Figure  23.  Further  suppose 
that  the  subgrade  has  a  modulus  of  23,000  psi  (approximately  CBR  15),  and  that 
the  design  aircraft  has  an  equivalent  single  wheel  load  of  5600  lb.  Entering 
Figure  23  with  a  subgrade  modulus  of  23,000  psi,  the  5600  lb.  line  is  inter¬ 
sected  at  100,000  coverages.  Hence,  the  proposed  pavement  should  sustain 
100,000  coverages  of  the  design  aircraft  before  developing  a  1  in.  rut. 

Further  full-scale  testing  to  include  instrumentation  of  the  sand  grids 
to  measure  the  lateral  stress  distribution  inside  the  cells  and  the  vertical 
stress  distribution  at  the  subgrade's  surface  will  allow  appropriate 
refinement  of  the  analytical  model. 

Extensive  exercising  of  the  design  model  will  allow  the  development  of 
complete  design  charts  which  can  be  used  directly  for  sand-grid  pavement 
design. 
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APPENDIX  A 
PROGRAM  LISTINGS 


-y~  -3  V'  MAE 


V®1 


w 

i 

w. 


A 

A 

i<8 

iL 

ffl 
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\V| 

S 


4 

?w’ 

m 

»Jya 

fta 

M 


& 

*s 


?o0-?Av  CAL  =rp  calcol/tin’g  tie  stest  or  "he  got  ape 

DIMENSION  R  W1-)  ,32  ( 1  2)  ,  G  1  (12)  ,0?  (12)  ,  W  (1)  , 3ZA  M  1  (12)  ,  3E  AM  2  (I7)  , 
r^iD  !  M2)  ,5  3102  (  T2) 

Q  *  1"*l  ft  *  \  #  *1  /  ▼  \  '  T  —  1  1  \ 


:i»  1  10)  {'•-SAMI  (I)  ,1=1/12) 

VPI"1  (3,  120)  i’3E  AM  2  (I1  ,1=1,12) 

SPITZ  P,  2  10)  O'5 IE  1  (I)  ,1=1,12 
V?JZ?l3t22 O]  (02102  T  ,1  =  1,  12 
A1=1.)5*1.75 
00  1  1=1,12 

R  1  (I)  =  3Z  1  (I)  *W  (1) 

«2  (I  =  BEAM 2  (I  *W  (2) 

0  1  (T  =  G P. T D  1  (I)  *1.0 

G?  r '  =g?id?  (i  *i.o 

CONTI  N7E 

W?.!"*1  (7,3  10)  (31  (I)  ,1=  1,  12) 

WRITE  (3  *  320  W32  I  .1=  1,  12[ 

S  1  =  2.  *  f .  55*  B 1  ( 1)  +  .fl26*$  1  (2)  ) 

S5  1  =  51/,\1 

5?=.  *»*2  1  (2)  +  . 45*B  1  (  1)  +.22*31(2)  +  .47*B2(1)+.  00  22  *B2  (2)  +  1 .  *G  1  (1) 
+.047*G  1(2) 

S52=52/A  1 

5  7  = .  ^ 22*^1  (2)  +  .  077*31  (3)  +.  27* B2  (1)  +.013*B2  (2) 

SS  3  =  S  3/A  1 

5  4=.  3*1*  1  (7)  +.  27*32  (1)  +  .  01  3*B2  (2) 
ssa=su/AT 

55  =  .  17*81  (3)  +  .  13*31  (4)  +  .27*3?  (1)  +  .013*B2(2) 

555  =  5  5/A  1 

56  =  .  42*31  (4)  +.?8*B1  (5)  +.  13*32  (1)  +  .006*B2  (2)  +  1 .  *G2  ( 1)  + .  047*0  2  (2) 

556  =  5  6/  14.9 

S9  =  2.*(.  25*31  (2)  +  .0  2  56*32  (2) +  .55*G1  (2)) 

5S9=S9/A  1 

5  10  = .  13  5*3  1  (2)  +.  064  *B  1  (3)  +.  01  9*B2  (2)  +.  0  066*B2  (3)  +  .  1  49*B2  (2)  + 

.  406*G2  (2  +.  141*02(3) 

55 10=  S 1  0/A  1 
c  i  -'5*h  1  Mw 


1  a  1  *G 


55 10=  S 1  0/A  1 

1  1=.  ”'5*3  1  (3)  +.0?5fi*B2  (3)  +.  1U9*B2  (2)  +  .  547*3  2  (3) 

551 1  =  51 1 / A  1 

E  1?  =  .  14  7*3  1  (3)  +.  106*B  1  (4)  +  .  015*32  (3)  +  .0  1  1  *32  (4)  +.  3  13*B?  (3)  + 
.  754*32  (4) 

5^ 12  =  5 1  2/A  1 

5  13= .  35*31  (4)  +.  31*B1  (5)  +.  25*B2  (4)  +.  13*32  (2)  +.02*B2  (3)  + 

.  73*G1]4)  +.  1  1  *G  1  5  +.  15*3  2  (2)  +  .  02*G2{3) 

5317=51 3/ 12.25 
5  16  =  2  .*  (-  1  1*3  2  (2)  +.  04*3  2  (3)  ) 

5^16=51  6/A  1 

5  17  =  .  15*32  (3)  +.  1 1  *P 2  (2)  +.04*B2  (3) 


vm 

9 

% 

wl 

"'?!-=  "  1  1 
”  1  3  =  .  O'?**?  /  3\  4-.  ns*r.-? 

5  5  13  =  S1VA1  ' 

*22=?.*  .  15*0?  (?) 

SS22=S?2/A  1 
C??=.  09*22  H)  + 

S523=  52  3/A 1 

524=.  a? *'■2  (4)  +.  33*3  2(5} 

♦  .23*S2]4 

5524=S24/12.25 

*29=2.*  (.  09*32  (  3)  + 

53  2  9=  S2  9/ A 1 
V7I”F  (3.  109)  s  1 . 551 ,  S2 , 552.  S'*,  553. 

♦  5 10,5S 10,511  ,SS  11 ,5  12 
"  12  ,551"-,  513 


(9)  +.  1  1*2  2  ('’)  +.  204*52  (?) 

96*32(4)  +.  15*3  2  {  3) 

.  12*32  (3)  +.  03*32  (4)  +  .36  *G  2  (3)  + 


06 *B2  (4) ) 


ion  pormatc// 


+ 

+ 

+ 

+ 

♦ 

+ 

+ 

+ 

+ 

+ 

+ 

♦ 

+ 

♦ 

♦ 

♦ 

♦ 


1  10 

120 
2  10 
220 
310 
320 


FOPflA"'  1M1 
POPMATf//, 
FORM  AT  (//, 
FOFM AT I  //, 
FORMAT  //, 
FORMAT  ’//, 

END' 


51 

52 

53 
S  4 
5  5 
56 
co 

5  10 
5  1  1 
5  12 
5  1  3 
5  16 
517 
^  1  8 
529 

523 
c  2  4 
5  29 
//  ' 


fW.3 
F  10.  3 
F  10.  3 
F  1  0 .  3 
*10.  3 
F10.  3 
no.  3 

F  1 0 . 3 
F  1  0  .  3 
3 


F  1  0 
F  1  0 
F  10. 

F  1 0 . 

F  1  0 . 

F  10. 

F  10 
F  10. 

--.710. o, 
BFAM1  LOA 


3FAM2  LOAD  = 


GRID  1 
on  D2 
BEAM  1 
BE  A  M2 


LOAD  = 
LOAD  = 
OF 
OF 


55  18 

3  X  , 

3Y, 

3X, 

3X, 

3X, 

3X, 

3X, 

3X, 

3  X , 

3X, 

3X, 

3X, 

3X, 

3X 

3  X , 

3X, 

3X, 

lh 


#  ^  2  2 
55  1 
S  5  2 
S  5  3 
55  4 
SSr> 

S  56 
5S9 
55  10 
55  1  1 
SS  12 
SS  13 
55  16 
SS  17 
S5  18 
SS22 
55  73 
SS  24 
SS  29 
12' 


t 

, 12[F2. 1 


54,554,55 
, SS 1 2, 5  1  ? 
,5522.52  3 
'  F^O.  3 
F  1  0.  3 
F  1  0.  3 
F  10.  3 
F  10.  3 
F  1  0.  3 
F  10.  3 
F  10.  3 
F  1  0 .  3 
F  10.  3 
F  1  0.  3 
F  1  0.  3 
F  10.  3 
F  1  0.  3 
F  1 0.  3 
F  1 0.  3 
F  1  0.  3, 
-  ,  F  1 0.  3) 


5S5. 56 ,SS6, 59,559, 

SS 13,S  16,5S16, 

5S23, S 24, 5524, S29, 5529 
3X,/, 

I 

/, 

I 

'A 

/. 


ij  -  ■  ,  i 

3.2  IS. 
6.4  IN. 


1 2  (F7 .  1 
”  LOAD 
LOAD 


11 


=  ',12  (F7.  11  1 
=  M2(F7.  ill 
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P^03r\M  5&R*t> 


n“Z5“  revision  -  hay  25  ,  1995  —  modified  by  chant, 

UNDER  "HZ  D  I?  ZCTON  OF  DP.  P.  LAMBS. 


??OGeAM  GSAND  IF  MCDI“IFD  RASED  ON  I.  !i.  SMITH  PROGRAM  5.0 

(PLANE  5"  FAIN  0-  SLAVIC  FOLIO  'JSJKc  4-NOOZ  DU  ADRTL  A^ZP.  AL  ELEMENTS) 

T':Z  MODIFICATIONS  ARE  FOLLOWS  : 

1.  ECHO  TN?'l"  P\~A 

2.  CAN  HANDLE  DIFFERENT  E,V  VALUES  *OR  EACH  ELEMENT  0°  LAYER 
9 .  CAN  HANDLE  PI FT” FEN"  ELE-MSM'"  SIZE 

4.  CAN  ASSEMBLE  GLOBAL  STIFFNESS  MATRIX  IN  N*N  MATRIX 

5.  CAN  HANDLE  ''MOWN  DISPLACEMENT  BOUND APY  CONDITION 
5.  CAN  CALCULATE  “HE  FORCE  (REACTION)  VECTOR 

7.  ?RIN“  S“FESC  AND  STRAIN  AT  SAME  LINE  THROUGH  NEW  SUBROUTINE 
PE INTW  “OP  EACH  GAUSSIAN  DOIN'"  OP.  PRIN“  LOAD  DEFOR.  FORCE 
AS  TABULAR  FORM  FOR  EACH  DEGREE  OF  FREEDOM. 

3 .  CAN  I“EPA“E  PSUEDO-NCNLINER  PPOBLEMBY  SEATING  rnHE 

CONVERGENCE  CRI“ERION  AS  1/10  THE  DIFFERENCE  OF  THE  SHEAR 
S“R  AIN  RA“IO  (G/GMAX)  BETWEEN  'PVO  I’T’SRAmIONS. 


REAL  Z,  V,  AA,33,K  I  ,  K  2,  D3T,  Ql’O? 

TN“EGE?  I, J,K,L,P,Q,NXE,NYE,H,N,W,CDMAX,3H,NL,MN, 

*10’’,  FOD,NODOF,R,~,GP,  IDEF , ICOOFD, I J AC, I JAC 1  ,IDERIV, 
*TDEE,IBEE,TDt|EE,T3T,IBTD3,IKK,INF,ISAHP 
REAL  DS“  (0, 3)  ,  SAMP  (7,2)  ,  COORD  (8,2)  ,  JAC  (8,2)  , 

*7  AC  1  ( 8,  P)  ,  DER  (8,4)  ,  DERI  V  (8, 4)  ,  BEE  (8,8)  ,DBEE  (8,8)  , 

*3TDP  (3,5)  ,KM(8,fl)  ,  ELD  (8)  ,VOL(8)  ,EPS  (2)  , SIGMA  (3)  ,BT(3,3)  ,FUN  (4) 
TN^EG"?  G[8) 

AL^E”  HEX'"  '"WO  CARD'5  '’'O  CHANGE  PP.0  3L EM  SIZE 

PEAL  BE  (  10  0  0)  ,  LOAD:,  (100)  ,  8IGE  (100,100)  ,FORCE(10C)  ,GDIS?(100)  , 

*  B“GXM  (  100,  100)  ,VAL  (IOC)  ,  RADIUS  (100)  ,  DEPTH  (10  0)  ,  PRO?  (10  0,2)  , 

*  GSTRN (23)  , PATIO (23)  ,3TEESV (22)  ,XO (22)  , G RATIO  (22,  10) , 

*  AV GS“N  (2  2)  ,  LOAD  (  10)  ,  GMAX  (22)  ,  LO ADD  (1  00)  ,  DISP  ( 100)  ,  FORCES  ( 1  00) 
INTEGER  NF  (  100,2)  ,  IBIGK  ,FIX  ,NO  (  100)  ,NLOAD  (1C) 

D  AT  A  IDE?, ICOOED, IJ AC, I JAC 1 ,TDERIV, IDEE/6* 8/ 

P A“A  I3EE,ID3EE,IB“,  IB'T'DB,IKM  ,DOF/6*8/ 

0  A“A  IS  AMP/7/, NODCF/2/,H/3/,T/2/ 

D  AT  A  GSTF.N/.  0000C  1 ,  .  0 00 0  1 , .  00  0  1  , .  000  1 5 ,  .  03  02,  .  0003  ,  .  000  4 , .  0 00  5 , 

*  .0006,. 00 07,. 0003  ,. 0  00 9 ,.00  10,. 00  15, 

*  .  00  20, .0030,. 00  40, .00 50,. 0060,  .0070, 

*  .0080, .0090, .0100/ 

DATA  RATIO/1. ,. 95, . 74,. 67, .^2,. 53, .43, .43, .39, .36, .34, .32, .30, .27, 

*  .  19,.  15,.  12,.  10  5, .09, .08, .074,. 065,. 0  6/ 


ALTER  NEXT  CARD  wO  CHANGE  PROBLEM  SIZE 
DATA  TNF/100/, I3IGK/100/, EPSLON/. 1/ 

INPUm  AND  INITIALISATION 
READ  (1, *)  N  X  S,  N  YE,N,W,NN,5N,NL, FIX, G° 

A-4 


(3,  100)  ‘JXE,*IY7,M,N,F'I,  ?N,NL  ,F"',G? 

100  7o*»m~  ( i»  1 ,/,  •  number  o”  fle me v""  im  x  -  nrrrio"  = 

♦  1  vUr'DSF  OF  ELEMENTS  IN  Y  -  DIRECTION  =  ',14,/, 

♦  *  "*0 "" \L  NUMBER  O?  77  EE'LOMS  I 'I  M"SH  =  ',14,/, 

♦  •  HAL-  3  AMD  WTD”'J  =  *,I4,/, 

♦  t  ^rrMBFP  07  NODEG  I«  MESH  =  ',14,/, 

♦  •  NUMBER  07  RESTRAINED  MODES  I'l  MESH  -  •  ,14,/, 

♦  •  NUMBER  OF  LOADED  FREEDOMS  =  ',14,/, 

+  .  0?  FIXING  PRESCRIBED  DISPLACEMENT  =  '  ,14,/, 

♦  •  NUMBER  OF  GAUSSIAN  IN^FPGR  AmION  ORDER  =  ',14,//) 


N0D=D07 /MODOF 
COM  AX=U ♦ 1 
R=N*CDM AX 
M  =  N  X  E  +  1 

R  SAD  ( 1 ,  *)  (RADIUS  (I)  ,  T  =  1  ,  M) 

WRI^S  (.1,  1  10)  (RADIUS  (I)  ,1=1, M) 

110  FORMA'"  (/,  •  RADII  07  "HE  MESH  LINES  =  »,2?F5.1) 

M=UYE+  1 

READ  {  1,  *)  ("E?mH  (I)  ,1=  1,  M) 

WRITE  (3,120)  (DEPmH  (I)  ,  1=  1 ,  M) 

120  FORMA'"  (/,  •  DS?mH  OF  THr  MESH  LIMES  =  ',1176.2  ) 

CALL  GAUSS  (SAM?, ISAM?, G?) 

CALL  FORMMF  (NF, IMF, NN , MODOF ,P N) 

ELEMEN m  STIFFNESS  INTEGRATION  AND  ASSEMBLY 

WRITE  (3,14  5) 

146  FORMA”  (  1 H  1 ) 

ITER  =  0 

1  I’ER  =  -rmp-j  +  y 
CALL  JT1LV7C  (BK,E) 

CALL  NULVEC  (LOADS,  N) 

CALx.  NULL  (BIGK,TPIGK,  M,  N) 

CALL  MUI.L  (3IGKM,  IPIGK,  N,  N) 

CALL  N'T  LI.  (DEE  ,  IDEE,  H,  H) 

DO  3  P»1,  NX  2 
DO  3  0=  1 , NYE 

NUM3ER  =  NYE  *  (P  -  1)  +  Q 
IF  (ITER  .EQ.  1)  THEN 

READ  {  1 ,  * )  STPFSV  (NUMBER)  ,  KO  (NUMBER)  ,V 
S  TR  ESC  =  ( 1 .  ♦  2 .  *K 0  ( MB  MB  SF ) )  /3  .  *  SmPES7  (NTJM  BSP) 

GMAX  (NUMBER)  =  5500.  *  (STRESO)  **  0.5 
E  =  2.  *  ( 1 .  ♦  V)  *  GMAX  (NUMBER) 

WRITE  (3,50)  ITER,  NUMBER, KO  (NUMBER)  ,  NUMBER ,  STREP  V  (  NU  MB  EP)  , 

*  •  STPSSO, GMAX  (NUMBER)  ,E 

rr  <?T7 

«  .J  ■  .  i 

ITER  1  =  ITER  -  1 

7-P7c:o  =  (1. +2.  *KO  (NUMBER)  ) /3.  *  STRESV  (NUMBER) 

GG  =  ”  500  .  *  GFATIO  (NUMBER,  ITER  1)  *  ( STRESO)  * *0 . 5 

E  =  2.  *  (1.  +V)  *  GG 

U  ^  I  t  7  (3,55)  I”  ER  ,  NUMBER,  GRA^IO  (NUM3EP,  ITER  1)  ,  STPES  V  (NU  M3ER)  ,* 

*  STFE  SO, GG , E 
ENDIF 

rO  FORMA” (/, »  TTE?= ,,I3,2X,,KO(,,I2,')  =  • , F9 . 2, 2X, • S V  ( • , 12 , 1 )  =  '  , 


A- 5 


mm 


w-/ 


MMPIU  »  'J  IU»PlVIUflJ  rj  ITJ  V|l|ri  f '»V  '  V  >  r«  "■>■•• 


v*’w^  *v"  ■  vh'Vt  v  ■  \r 


—  o  7  ^  y  »  7  ^  r  »  c  n  *"> 

ro-.:u"  (/,  •  T.-Tr=  •  ,  :3,~:',  'g?a" 
2Y, • S0  = '  ,F?  . 2 ,2  X,  ' 
CALL  NH L L  (KN, IKS, DO", DO") 

CALL  FOP  HD  (DFE, IDEE, E,  7) 

TALL  V  G  E  C  N  Q,  N  X7  ,  S  YF  ,  CCO  HO  , 

DO  4  1=  1  ,G? 


1=  '  ,^3 .2 ,2? ,  '  ?=•  ,~i.2) 

10  ( ’  ,1?  ,  •  )  =  •  ,  F9 . 2  ,?  <,  '  SV=  •  ,  F9  . 
3  = ' ,F9. 2, 2X , • E  = ' , F9.2) 


IC003D, G, NF,  INF,  F  A  DINS,  DEPTH) 


D^  4  .7=1,3? 

X  1  =  [I,  2) 

K2  =  SAX?  (J,  2) 

CALL  FOPNI. N  (DF? ,  I  DF?  ,  Fll  N  ,  S AH?  ,  I S ANP , I  ,  J) 

CALL  N»"MnL (DS?,IDSR, COOFD, ICOOPD , J AC, I JAC, T, NOD, T) 
CALL  "W0BY2  (J  AC,  I  JAC,  JAC  1  ,  IJ  AC  1 ,  DE") 

CALL  HAmHrTL  (G  AC  1,  TJ  AC  1,  DSP  ,  IDEP  ,  DSP  IV  ,  IDE  PI  V,  T,  7,  NOD) 
CALL  N'JLL  (BEE, LESS, H, DOF) 

CALL  r0R  N  3  (3r>S,I3EE,DFPIV,TDF?.IV,VOL,NOD) 

CALL  H  ATHOL  {DEE  ,  IDF  E,  BE  Z,  13  EE ,  D3EE,  ID  BEE,  H,  H,  DOF) 

C  ALL  HA"? AN ( 3T , I3T, BEE, IB  EE , H, DOF) 

CALL  HA"HCL (3" , IB" , DBEE , I DBEE , B^DB, TB"DB , DOF , H , DOF) 
Q'J0"  =  DET*K1*K2 
DO  5  K=  1  , DOF 
DO  5  L  = 1 , DOF 


2, 


S  B"  D3  (K,L)=3"D3  (K,L)  *QDO'T' 

CALL  HATADD  (K  H  ,  IKN  ,  B^DB  ,  IBTDF. ,  DOF,  DOF) 
4  CONTINUE 

CALL  ^HBIGK (3TGK, I3IGK, KH, IKH,G, DOF) 

CALL  FH3IGK ( BIGKK, IBIGK, KM , IKM, 3, DOF) 

CALL  FORHKV  (BK,  K.H,  IK  K  ,  G  ,  N,  DO  F) 

1  CONTI N  IS 


EQUATION  SOLUTION 


IF  (ITEP  .  FQ.  1)  THEN 
DO  i  1  1=1, "IX 

P FAD  {  1 ,  *)  NO  {I)  ,  VAL  (I) 

11  WP.I"*  (3,  170)  NO  (I)  ,  VAL  (I) 

ENDIF 

170  FO?H  A"  (/,  •  ?E  ESCF.IB  ED  DEFOP.MA"ION  A^  «  OF  DOF  (*,12,')  =  * ,  F  T  2.  4 ) 
DO  12  1=1, FIX 

3K  (NO  (T)  )  =  BK  (NO  (I)  )  ♦  1.0E12 

12  3IGKN  (NO  (I)  ,NO  (I)  )  =  BIGK H  (NO  (I)  ,  NO  (I)  )  ♦  1.0S12 
CALL  BANPF.D  (3K,N,CDNAX) 

DO  1?  1=1, FIX 

13  LOADS  (NO  (I)  )  =  BK  (NO  (I)  )  *  VAL  (I) 

I*  (I" S3.  .  SQ.  1)  "HEN 

DO  6  1=1, NL 

BEAD  (  1,*)  NLOAD  (I)  ,LOAD(T) 

WFI"E (3,  18H)  NLOAD  (I)  ,  LOAD  (I)_ 

8  LOADS  (NLCAD(I)  )  =  LOAD  (I) 

END  IF 

180  FO°H  AT  (/,  •  PPSSC?  I3ED  LOAD  AT  *  OF  DOF  (*,T2, ’)  =  *  /F 1 2^- 4) 

I ?  ( I T EF  .  GE.  2}  "HFN 
DC  6  8  1=1,  NL 

66  LOADS  (NLOAD  (I)  )  =  LOAD  (I) 


DO  2  10  1=1,0 

LO  ADD  D  =  LOADS  (I) 

■)  rcMTI?!'l" 

CALL  BACK  1  px,  LOADR,  ?J,CD.M  AX) 

C  *•  L  L  9ACX2  (BK,  LOADS,  0,  COM  AX) 

00  0  1C  1=1,  'I 

D  ~  3?  (I)  =  LOADS  { I ) 

OK  c^'PI  on  0 

C5LL  MV  MOL"  C3!  OX,  19  IS  X,  LOADS,  N,  0,  FORCE) 

00  4  10  1=  1 ,  N 

FORCES  (I)  =  FOPCE(I) 

4  10  CONTI  N’J  E 

CALL  XVMni.T  (3IGXM, I3IGX,  LO  ADS  ,  N,  N ,  FORCE) 

W?ITE  (  1,550) 

550  FOP ‘I  A" (  1H  1 ,  *  NF  KBE*3  *  ,  7X  ,  '  LO  AD  *  ,  K  X,  *DIS?LACEMSNm,,6X,  'FORCE*  ,/) 
DO  5  0  0  1=1,?! 

WT,ImE  (3,  ■’0  0)  I,LOADD  (I)  ,DISP  (I)  ,  FORCES  (I) 

500  ccnti:dje 

700  FORMA7' (15, 5X,  7  {E12.«,  5X)  ) 

PEAC"  =  0. 

DO  900  1=1, FIX 

REACT  =  PE  AC7,  *  FORCES  (MO  (I)  ) 

900  CONTINUE 

W^ITE  (3,850)  ITER , REACT 

850  cO  P.M  A  "  ( /,  *  REACTION  FO^CE  FOP  I^ER.  (’  ,12,  *)  =*,F10.2,/) 

RECOVER  ELEMEN"  STRAINS  AND  STRESSES 
AT  ALL  GADSSIAN  INTEGRATION  POINTS 


\ Y y  =  ?.’ yf  *  NYE 
DO  272  1=1, NXY 
*  7GSTN  (I)  =  0. 

GDIS  ?  (I)  =  0. 

222  CONTI  N"IE 
DTrMAX=0. 

DO  7  ?  =  1 ,  N X S 
DO  7  q= 1 , NYE 

CALL  NNLL  (DEE, IDES, H,  H) 

C'LL  (DEE,  IDES,  E,  V) 

CALL  VS EON (n, Q , NXE , NYE , COO RD, ICOORD, G, NF, INF, PA DIOS, DEPTH) 
SDN  =  0. 

DC  8  1=1, G® 

DC  8  J=  1  ,GP 

CALL  FOR  MLN  (DER,  IDER  ,  FU  N,  S  AM?  ,  ISAMP,  I ,  J) 

CALL  MA"MUL  (D  ER  ,  IDER  ,  COORD ,  ICOORD  ,  J  AC  ,  I J  AC,  T,  NOD ,  T) 

CALL  "1*0372  (JAC,  IJAC,JAC1,TJACJ,DE’T') 

CALL  MATMUL  (JAC  1 , 1  JAC  1,  DEP,  ID  EP,  DER  IV  ,  IDEP.IV  ,  T ,  7,  NCD) 
CALL  NULL  (BEE,  TREE,  H,  DOF) 

CALL  EO^MR (DEE, TREE, DER IV , IDSRI V , VOL, NOD) 

DO  9  M= 1 , DOF 

IE  (G  (1)  .  EQ.  0)  ELD  (M)  =0. 0 
0  IF  (G  (“)  .  NE.  0)  ELD  (M)  =  LOADS  (G  (M)  ) 

CALL  MVMUL"  (DEE, I  BEE,  ELD,  FI,  DOF, EPS) 
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wrrwr*w^mrwi 


f 


CALI  M7MUL"  (:■”,  IDEE,  E’%  H,  :!,  SIGMA) 

SUM  =  SUM  +  *??>  (3) 

number  =  ny"  *  {?  -  1}  ♦  q 

3  CONTINUE 

AVG^"N  (NUMBER)  =  SUM  /  4. 

GD IS?  (MEMBER)  =  SUN  /  4. 

DO  3  1  K=  1,2.3 

I*(GSTDN(K)  -  ABS  (A VIST N  (NUMBER)  )  )  31,82,32 
81  CONTINUE 

33  K 1 1  =  K  -  1 

I~(K11  .  EQ.  0)  GO  "0  3  3 

SLO?E=  (PATIO  (K)  -RATIO  (K1  1)  )  /  ALOG  (GSTRN  (K)  /GSTPN  (K  1  1)  ) 

=  SLO^E  *  ALOG  (A3S  (AVGSTN  (NUMBER)  )  /  G3TRN  (K  1  1)  ) 

G^rm  (NUMBER,  rE?)  =PA'T'IO(Kl  1)  -ABS  (mT) 
a3  IF  (K  1  1  .EQ.  0)  GP  ATIO  (NUMB  EP  ,  ITER)  =RATIO  (1) 

IP  (ITS?  .GP.  2)  "HEN 
I"ER 1  =ITE?  -  1 

TIFF  =  A3S(S*>  ATIO  (NUMBER,  I"ER)  -  GR AT  10  ( NUN  BER ,  IT ER  1)  )  / 

*  ABS  (GI?  A"IO  (NUM3ER,  I"ER  1)  ) 

IF  (APS  (DIET)  .  GE.  DTFMAX)  DIFMAX  =  APS  ()  IFF) 

END  IF 

MPT'" E  (3  ,1150)  NUN3ER,  AVGSTN  (NUMBER)  ,  NIT MB EP.  ,  G RATIO 

*  (NUMBER, ITER) ,DIFF, DIFMAX 

1150  FORM  A1"  (  1  **  ELEMFN"  =  ‘,I3,2X,’  AVG.  SHEAR  SmRAIN  =»,F16.8,2X, 

*  'GPATIO  (*  ,12,  »)  =•  ,F6. 2,2X,*  DIFF  =», F 1 0. 8, 2X, ’ DIFMAX  =‘,?16.8,/) 
3  CONTINUE 

IF(I"E?  .EQ.  1  .OR.  DIFMAX  .GT.  SPSLON)  GO  TO  1111 
0339  WRITE  (3,  1250)  I*"ER 

*0  FORMA"  (/,  •  ITERATION  =  ',’‘3) 

STO? 

END 

rMB?00"  INE  GEOMF"  (P,Q,NXE,MYE,  A  A,  BB ,  COORD,  ICOOR^,  G,  NF,  I NF) 

NODAL  COORDINATES  AND  STEERING  VECTOR 
FOR  A  RECTANGULAR  MESH  OF  4-NODE 
QUADPI L  AmER AL  PLANE  ELEMENTS 
NUMBERING  IN  THE  Y-DIRECTION 


IN" EGER  G  { 1)  ,  INF,  *’F  (INF,  1)  ,  P,  Q ,  NX E,  NYE ,  ICOORD ,  AO,  AL,  AM,  AN 
PEAL  COORD  (ICOORD,  1 )  , AA , BB 
AO=  (P-1)  *  (NYE*  1 )  +Q 


AL=  AO* 1 

\H=r>*  (NYE+  1)  +  Q 
AN=  AM  +  1 
G  ( 1 )  =  NF  (AL,  1) 

G  (2)  =NF  (AL,  2) 

G  (3)  =  NF  (AO,  1) 

G  (4)  =  NF  (  AO,  2) 

G  (5)  =NF  (AM,  1) 


G  (6)  =*?  (AM,  2) 

G  (7)  =  NF  (AN,  1) 

G  (3)  =  NF  (AN,  2) 

COO  PD  (1,1)=  (P-1)  *  An 
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crop?  (1,2)=  jmye-o)  *rs 
cooao  p,  i)  =  (?-i)  * \ a 
COORD  (2  ,  '*)  =  (NYE-Q+  1}  *F3 
COO^O  (2,  1)  =?*  A  A 
COCR?  (3,2)  =  (NYE-Q  +  1)  *DB 
COO  PC  (4  ,  1)  =*>*  A  A 
C^ORP  (4,2)  =  (NYE-Q)  * 3 P 
P  E'"H?  N 
E'JD 

S'JPPO’J-TMZ  N’lLVEC  (VZC,  V.) 

PEAL  VEC  ( 1) 

IN-EGER  I, 'I 
on  1  r=i,'j 
1  VZC  (I) =. 0 
^Z'UEN 
EMD 

SUPPOnrE  'HILL  (A,  I A  ,  M,  N) 

IN-EGRR  IA,M,N,I,J 
PrAL  A  (IA,  1) 

DO  1  1=1,  M 
DO  1  J=  1 ,  *} 

1  A  (I,  J)  =0.0 
RE-U2N 
END 

SUB  RO'J- 1 NZ  POP ’ID  (DEE,  IDEE,  E,V) 

STRESS  STRAIN  MATRIX  FOP  PLANE  ELASTIC  STRAIN 

IN-EGZR  IDEE, I, J 

PEAL  DEE  (IDEE,  1)  ,E,V,  VI  ,VV 

V1=V/(1.0-V) 

W=  {  1 ,0-2-0*V)  *.  5/(  1 . 0-V) 

DEE  (1,1)  =  1.0 
DEE  (2 ,2) =1 . 0 
DEE (3,3) =VV 
DEE  ( 1 ,2) =V  1 
DEE  (2,  1)  =V  1 
DO  1  1=1,3 
DO  1  J  =  1,3 

1  DEE  (I,. 7)  =DEE  (I,  J)  *E/(2.  0*  ( 1  .0  +V)  * VV) 

PE'-tJPN 

END 

SUBROU'?  IMP  GAGES  ( S  A  MP  ,  I S  A M P  ,  G  P) 

GAUSSIAN  QUAD? AmUP  E  ABSCISSAE  AND  WEIGHTS 

7  N"  EGEP,  ISAMP,GP 
PEAL  SAMP  (ISAMP,  1) 

GO  -0(1,?,  3, 4,  S,  6, 7)  ,GP 

2  SA^P  (1,  1)  =  1  -  /SORT  (3.) 

SAM?  (2,  1)  =- *7 AMP  ( 1  ,  1) 

SAMP  ( 1 , 2)  =  1 . 

SAM?  (2,  2)  =1. 

GO  TO  1 


->  c\-‘?  (1,1)=.  "’♦so?"  (15.) 

SA  v?  (2,  1)  -.0 

p,  i)  =-sam*mi,  i) 

SA'*?  '  1,2) 

SAM?  (2,  2)  =8./*?. 

SAMP  (3,  2)  =  5A.MD(1,2) 

GO  -0  1 

u  SAMP  (1,  1)  =  .R6  1136311594053 
8  A"?  (2,  1}=.  3?  99  81043584  856 

samp  (3,  i)  =-sam°  [2,  i) 

S\M?  (4,  1)  =-~AM°  (1,  1) 

SAwr>  (1,2)  =•  3478 54R4  5  13-454 
GAMP  {2,2.)  =.  6521451548  62  546 
SAMP  {3,  2)  =  S AM*  (2,7) 
sAM?  (4,  7)  =SAM?  { 1  /  7 ) 

GO  "0  1 

5  *■  AM?  ( 1 ,  1)=.  906  179845938664 
S  AMP  (2,  l)  =.938469310105683 
SAMP  (3/  1)  =.  0 

SAMP  (4,  1)  =-SA.M?  (2,  1) 

SAW?  (5,  1)  =  - S  AM?  { 1  /  1 ) 

SAM? (1,2)=.  23692688  50  55  189 
s  AMP  ( 2,  2)  =.  4  7  86  2867  04  99  366 
SAMP  (3, 2) =.  56888388  33  88889 
SAMP  (4,  2)  =  S\M?  (2,  7) 

SAMP  (r>,2)  =  SAMP  (1,2) 

GO  '”0  1 

6  SAMP  (1,  1)=. °32469514203152 
SAMP  (2,  1)=.  66  1209336466265 
SAMP  (3,  1)  =.  2386  19186083  197 
SAMP  (4,  1)  =- S  AM*  (3,  1) 

S  AM?  (5,  1)  =-  SAM?  (2,1) 

SAM?  (6,  1)  =-  SAM?  ( 1 ,  1) 

SAM?  (1,2)=-  17  1324492379170 
SAM?  (2,  2)  =.  369’76i5"T3043  139 
SAMP  (3,2) =. 46"9  1393  457269  1 
SAM?  (4,  7)  =  S  AMP  (3,2) 

SAMP  (9, 2)  =  SAM?  (2,2) 

S  A’lP  (6,  2)  =S  AM?  ( 1 ,  .2) 

GO  TO  1 

7  SAMP(1,  11  =.  949107912342759 
SAMD (2,  1)  =  .74  1531  16  5599  394 
S  AM?  (3,  1)=.  405845151377397 
SAM?  (4,  1)  =.  0 

SAM?  (5,  1)  =- S AMP  (3,  1) 

SAMP  (6,  1)  =-SAMP  (2,  1) 

SAMP  ( 7 ,  1)  =-SAM?(1,  1) 

SAM? (1,2)=. 129484966 168870 

s  AM5  ( 2,  2)  =  .  2797053  9  1489  277 
5  AM* (3,  2) =.  38  18300  50  505  1  19 
SAMP  (4,  2)  =.  4179  5918367346  9 
SAMP  (5,2)=SAM*  (3,2) 

SAM?  (6 ,  2)  =? AMP  (2,  ?) 
SAM?(7,2)  =  S  AM*  ( 1 , 2) 
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i  Cr*i"‘rT" 

•>  T“*»|  ”3  >1 

END 

ST»  ROUTINE  70  = 'INF  (HF,  INF,  NN, NODOF,  RN) 

NODE  FF  EEDOM  ARRAY  FOP  NOPE 
TH  AN  ONE  FREEDOM  PEP.  NODE 

IV-FGE?  INF ,  NF  (TKF ,  1)  ,  NN,  NODOF,  FN,  I ,  J,  K  , L 
DO  1  1=1, NN 
DO  1  J=  1 , NODOr 
1  NF  (I ,  J)  =  1 

DO  2  1=1, FN 
?  F  AD  ( 1 ,  *)  K 
DO  3  J=  1 , NODOF 
PEAD  ( 1  ,*)  L 

3  IF  ( I- .  3Q.  1)  NF  (K,  J)  =0 
?  CONTINUE 

K=  1 

DO  4  1=  1 ,  NN 
DO  4  J=  1, NODOF 
IF  (NF  (I ,  J)  )  5,4,5 
5  NF  (I,  J)  *K 
K=K  +  1 

4  CONTINUE 
RETURN 
END 

SUBROUTINE  ’’ORMLN  (DSR,  IDER,  FfTN,  SAMP, IS AMP, I ,  J) 


LOCAL  COORDINATE  SHAPE  FUNCTIONS  AND  mHEI? 
DERIVATIVES  FOR  A  4-NODE  QUADRILATERAL 


INTEGER  IDER, ISAM?, I, J 

°Z  AL  OS’5  (IDER,  1)  ( i)  ,  R  AM?  ( IS  AM?  ,  1) 

REAL  E".A,XI,ETAM,E?AP,XIH,XI? 

E"'A=S  AM?  (I,  1) 

XI  =  S AMP  (J,  1) 

E'rAM=.?5*  (l.D-E^A) 

ETA  ?  =  . 25* (1.0+ZTA) 

XTM=.25*  (1.0-XT) 

XT?=.  ?5*  ( 1. 0+XI) 

FUN  (1)  =4. 0*XIM*ETAM 
rr7N  (2)  =4. 0*XIM*EmA? 

FUN (3) =4.0*XIP*E"A? 

FTTN  (4)  =4. 0*XIP*E?AM 
DEP  (1,1)  =-E"AM 
DE° (1,2) =-ETAP 
DEP  (1,3)  =  E'M? 

DEP ( 1,4) =ETAM 
Dp?  (2,  1)  =-XTM 
DEP  (2,2) =XI M 
DE?  (2,3)  =XIP 
DEP  (2,4) =-XIP 
RETURN 
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END 

SNBrn'r  t»|^  v\mMrL  f  A,  ~  A,  0,  :3,C,IC,L,M,N) 

7"r EOT?  TA,I3,Tr,L,M,N,T, J,K 

"**L  MIA,  1)  ,3(13,  1}  ,C(IC,  1)  ,X 
on  1  1=1,1 
00  1  J=1,N 
x=o.  0 

on  2  x=i,m 

2  X=X  +  \  ( I ,  K )  *3  (K,  J) 

C(I,J)='C 

1  ro STINTS 

T}  P'T’P  t>  'J 

END 

SyS’'OU'“TNE  "’W03Y2  (JAC,  1 .7 AC,  JAC  1  ,T JAC 1 , DE^) 

INTEGER  I JAC, I JAC 1 , K, L 

REAL  JAC  (I  JAC,  1)  ,  JAC  1  (T  JAC  1,1)  ,  DET 

OE^JAC  (1 ,  1)  *  JAC  (7,  2)  -JAC  (1,2)  *JAC  (2,1) 

JAC  1  (  1,  1)  =  JAC  (2,2) 

J  AC  1  (1,2)  =-JAC  (1,2) 

JAC  1  (2,  1)  =- J  AC  (2,  1) 

JAC  1  (2,  2)  =  J  AC  (1,  1) 

00  1  K=  1 , 2 
DO  1  L=  1 , 2 

1  JAC  1  (K,  L)  =  J  AC  1  (K,  I.)  /DE'r 
RETURN 
END 

SUBROUTINE  EORMB  (BEE,  IBEE,  DER  IV  ,  IDERIV  ,  VOL,  NOD) 

STRAIN- DISPLACEMENT  MATRIX  FOR  PLANE  STRAIN  (STRESS) 


t^ECE?  IBEE,IDERIV,NOD,K,L,M 

REAL  BEE(TBEE,  1)  ,DEPIV  (TOEF.IV,  1)  ,V0L(1) 

00  1  3=1,  NOD 
K=2*M 
L  =  K-  1 

VOL  (L)  =DE? I V  ( 1 ,  M) 

SEE  (  1,1)  =OERIV  (1,3) 

BEE  ( 3  ,  K)  =DERIV  (1,3) 

V0L(Y)=DERIV(2,M) 

SEE  (2  ,  K)  =DE^IV  (2,  M) 

SEE  (3,L)  =OEEIV  (2,M) 

1  CONTINUE 
RETURN 
END 

SUBROUTINE  FMBIGK  (BIOK,  IBIGK,  KM , IKM, G,  DOF) 


INTEGER  T3IGK,IK3,G(1)  ,  DOF ,  I,  J 
REAL  BIGK  ( I  3 1 G  K  ,  1)  ,  KM  {I  KM  ,  1) 

00  1  1=1, DOF 
IF  (G  (I)  .  EQ.  0)  GO  -0  1 
2  00  3  J= 1 , DOF 

It‘  (G  (J)  .  EQ.  0)  GO  m0  3 
4  ?TGK  (G  (I)  , G  (J)  )  =BIGK  (G  (I)  ,G  (J)  )  +  KM  (I,J) 


3  CONTI  TT- 
i  cc"!-: vm 

D  ?  U  *  V 

END 

SUB?r *T"“:iE  FORM KV  (  =  K  ,  KM  ,  IKK,  0  ,  N ,  DOF) 


ASSEMBLES  ELEMENT  MATRICES  INTO  SYMMETRICAL  BAND 
01, 03 AL  MATPIX  (SmORED  AS  A  VECTOR) 


INTEGER  0(1)  ,TKM,N,DO?,I,  J,CD,VAL 
PEAL  BK  (1)  ,  KM  (TKM,  1) 

DO  11=  1,D0F 

IF  (G  (I)  .EQ.  0)  GO  TO  1 

2  DO  5  J=  1/DOF 

I'®'  (G  (J)  .FQ.  0)  GO  TO  5 

3  CD=G  (J)  -G  (I)  ♦  1 
IF  (CD-I)  %4,4 

4  V  AL  =  N*  (CD-  1)  +0  (I) 

BK(VAL)  =BK  (VAL)  +KM  (I,  J) 

5  CONTINUE 
1  CONTINUE 

RETURN 

END 

SUB POUT  INS  MAmP  AN  ( A , I A , B, IB , M , N) 
INTEGER  IA,TB,M,U,I,J 
PEAL  A  ( I A  ,  1)  ,U  (IB,  1) 

DO  1  1=1,  M 
DO  1  J= 1 , N 
1  A  (J,  I)  =B  (I,  J) 

PFTURN 

END 

SUBROUTINE  M A” ADD  (A, I A, 3, T3,M ,N) 
INTEGER  IA,IB,M,M,I,J 
REAL  A  ( I A  ,  1)  ,3(13,  1) 

DO  1  1=1,  M 
DO  1  J=  1 ,  N 

1  A  (T  ,  J)  =A  (I,  J)  +B  (I,J) 

RETURN 

•  2ND 

SUBROUTINE  3ANPED  (BK,  L,  KB) 

REAL  BK  (1)  ,SUM 

INTEGER  L,KB,I,IL1,K3L,IJ,J,NK3,NI,NJ 

DO  1  1=2,  L 

TL  1  =  I-  1 

K3L=TL 1 +K3 

IF  (KBL-L)  3,3,2 

2  KBL=L 

3  DO  1  J=I,KPL 
IJ=  (J-T)  *L+- 1 
SUM=BK(IJ) 

M*B= J-K3+  1 
IF(NKB)  4,4,5 

4  NK3=  1 

5  T*  (NK3-II.  1)  6,  6,  8 


vV>  V.A 


*  W  “  T  - 


E  do  7  :;=';Kr>,TL? 

VI  =  (I  -N')'*L+  'I 

:;.7=  (J-M)  * l+  n 

7  <;rJM=  =1*1^- rjsr  'MT)  *BK  J)  /3K  (N) 

0  3X  {I  J)  =SU?» 

1  CONTINUE 

END 

c'TB?Otjn  INI  B^ok  1  (3K,?,7.,K3) 

?Z\L  B  X  (  1 )  ,  0  ( 1 )  ,  S  U  M 
TN^ZCZr  L ,  K  3, I, T1,NK3,JN,K 

p  (l)  =b(D/b:<(1) 

BO  1  1=  ?,  L 
SMM  =  R  (I) 

I  1  =  1-  1 
NKB=T-K3+ 1 
1 7  ( N  '<  3)  2,2,  1 

2  >!XB  =  1 

7  DO  4  K=  MK B , 1 1 
JN=  (I-K)  *L +  :< 

7UM=7UM-BK  (JN)  *R  (K) 

4  CONTINUE 

n  (T)  =  STT M/BK  (I) 

1  COM"  I NTJ  7 
RETURN 
END 

SU9P.0UTINS  PACK2  (BK,R,L,KB) 

R*‘»L  BK  (1)  ,R(1)  ,  STTm 
INTEGER  L,KB,JJ,T,I1,NKB,JN,K 
DO  1  JJ=2,L 
T=L- JJ+ 1 

p;tm=o.o 

I  1  =  1+  1 
***■9=1-1  +KB 
T9  (NXB-L)  7,3,2 

2  N'<3=  L 

3  DO  4  ;<=I1,NKB 
JM=  (X-I)  *L+ 1 

4  50M=SrJ;4*BK(  JM)  *B  (K) 

P  (I)  =  R  (!)  -SON/BE  (-) 

1  COM"!  NU  E 
RETURN 
END 

SUBROUTINE  .TUMULT  (B#  IN,  V ,  K,  L,  Y) 
I  V'ZGE?  IN  ,  K  ,  L ,  I ,  J 
P^L  «(IN,1)  ,  V  ( 1)  ,  Y  ( 1 )  ,X 
DO  1  1=1,  K 
X  =  0. 0 

DO  2  J=  1 , L 

2  X  =  X  +  N  (I  ,  J)  (J) 

Y(I)  =X 

1  CONTINUE 

END 
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~ 'i3 ro :ri': 7  ??rv-v  (v#l) 

REAL  V  (L) 

WRITE  (3  ,  1)  V 
1  ^ORMA'-f  It:  #  10^12.  4) 

r»;n 

S  UR  ROUTINE  PRINTI  (V,L) 

INTEGER  L 
->E\L  V  (L) 

DO  10  1=1, L 
W«ITT  (3,  1)  I,V  (I) 

1  FORMAT  (  1 1I  ,  •  NUKREP.  OF  FREEDOM  (•,12,*)  =  «,E12.4) 
10  CONTINUE 
RHTUPN 
HMD 


SUP  ROUTINE  VGHOM  {?,  Q,  NX H, HYF, COOPD, ICOORD, G, NF, INF, 

♦RADIUS,  DEPTH) 

NODAL  COORDINATES  AND  STEERING  VECTOR 
FOR  A  VARIABLE  RECTANGULAR  MESH  OF  4-NODS 
QUADRILATERAL  PLANE  ELEMENTS 


IMTEGEP  ?,  Q  ,  NX  R ,  NYE  ,  ICOOPD,  IN F,  G  ( 1)  , NF  (INF,  1) # AO, AL, AM, AN 
PEAL  COORD  (ICOORD,  1)  ,  RADIUS  (1)  , DEPTH  ( 1) 

AO= (”-1) *  ( N  Y  H  +  1 ) +Q 
AL=  AO-*- 1 

AM=  P*  (NYE  +  1)  +  Q 

AN=AM+ 1 

0(1)  =NF  (AL,  1) 

G(2)  =NF  (AL,  2) 
a  (3)  =NF  (AO,  1) 

G  (4)  =NF  (AO,  2) 

G(5)=N?(AM,  1) 

G(f)=NF  (AM,  2) 

0  (7)  =NF  (AN,  1) 

G  (B)  =NF  (AN,  2) 

COORD  { 1  ,  1)  =RADIUS  (P) 

COOPD  (2, 1)  =R  ADI  US  (P) 

COORD  (1,2)  =  DEP’rH  (Q  +  1) 

COORD  (4,2)  =DSPTH  (Q+  1) 

COORD  (2,2)  =  DEP^H  (Q) 

COOPD  (3,2)  =DE?TH  (Q) 

COORD  (3,1)  *  RADI  US  (P  +1) 

COORD  (4,  1)  RADIUS  (P+1) 

PE”UPN 

END 

**  ************************  ********  4>*  ******  **********  ******41*41*  ****** 
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APPENDIX  B 


SUMMARY  TABLE  OF  WEIGHTING  FACTORS 


m 


»ft 

•ft 

•ft 

*.•!< 


Square 


*r*'Tor.L<.- 


1  v-»>>3.)vV 


Corner 


Weighting  _  (Width) (Factor) 
Factor _ Area _ 

0.55 

0.55 

0.026 

0.026 

0.30 

0.45 

0.022 

0.047 

0.0022 

0.222 

0.077 

0.27 

0.013 

0.30 

0.27 

0.013 

0.17 

0.13 

0.27 

0.013 

0.42 
0.38 
0. 13 
0.006 
1.0 
0.047 


Contribution  neglected 

Same  as  Square  2 

0.25 

0.0256 

0.25 

0.256 

0.547 

0.547 


0. 185 
0.064 
0.019 
0.0066 
0. 149 


(Continued) 


(Sheet  1  of  4) 


r'.v'v:/' 


Weighting  _  (Width) (Factor) 
Factor  Area 


0.406 

0.141 

0.25 
0.0256 
0. 149 
0.547 

0.142 

0.106 

0.015 

0.011 

0.149 

0.313 

0.234 


0.35 

0.31 

0.25 

0.22 

0.13 

0.02 

0.77 

0.11 

0.15 

0.02 


Same  as  Square  3 

Same  as  Square  9 

0.11 

0.04 

0.11 

0.04 

0.15 

0.11 

0.04 

0.09 

0.06 

0.11 

0.04 

Same  as  Square  4 
Same  as  Square  1 1 

(Continued) 
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Square 

Corner 

Weighting  (Width) (Factor) 

Factor  Area 

21 

Same  as  Square  17 

'  22 

3,3 

0.15 

13,3 

0.15 

23 

3,3 

0.09 

3,4 

0.06 

13,3 

0.15 

> 

24 

3,4 

0.42 

3,5 

0.38 

13,3 

0.12 

13,4 

0.03 

14,3 

0.86 

14,4 

0.23 

25 

Same  as  Square  5 

26 

■ 

Same  as  Square  12 

27 

Same  as  Square  18 

Same  as  Square  23 


Same  as  Squa.c  6 
Same  as  Square  13 
Same  as  Square  24 


1st  sub index  of  corner: 

1  -  Beam  1 

11  -  Beam  1  from  the 

other  direction 

2  -  Beam  2 

12  -  Beam  2  from  the 

other  direction 

3  -  Grid  1 

13  -  Grid  1  from  the 

other  direction 

4  -  Grid  2 

14  -  Grid  2  f’-om  the 

other  direction 


(Continued) 
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Square 


Corner 


Weighting 

Factor 


(Width) (Factor) 
Area 


2nd  sub index  of  corner: 

The  corner  number  on  the 
beam  referenced  by  the 
first  subindex 


APPENDIX  C 


SUMMARY  OF  INPUT  AND  OUTPUT  FOR  SGRID  ANALYSIS 


K,  =  0.38 
<r,  =  54.5 

E  =  83010 

0.48 

14.45 

45828 

0.48 

14.45 

45628 

2.44 

1.10 

21805 

2.44 

1.10 

21805 

mm 

BfjB 

3.50 

0.27 

12717 

P|Ij 

■ 

K(  --  Q40 

0.52 

0.52 

1.94 

1.94 

3.50 

3.50 

3.50 

3.50 

3.50 

3.50 

Oy  ~  33.S 

13.35 

13.35 

1.65 

1.65 

0.82 

0.82 

0.82 

0.82 

0.82 

0.82 

E  =  60577 

44743 

44743 

24329 

24329 

22026 

22026 

22026 

22026 

22026 

22026 

I  I 

0.025"  0023 


0.0 


I 

0.011" 


?  ?  ?  f  »  » 

0.005"  0.004"  0.003"  0.003"  0.0025"  0.002" 


E  =  51152 

18896 

16172 

15787 

7517 

8331 

9486 

10549 

12196 

11582 

10885 

18350 

17776 

19522 

15266 

16015 

17375 

19236 

21119 

18492 

t 

33.5lb 


i  i  A  ......  • 

I  zJlb  8.5ib  ,-7'b  -0.9lb  -0.2lb  -0.2lb  -0.2lb  -0.8lb  +1.2lb 


60lb 


33lh 


CASE  2-1  BEAM  1 


C-2 


.vv 


3.50 

0.27 

12717 

3.50 

0.27 

12717 

3.50 

0.27 

12717 

3.50 

0.27 

12717 

2.11 

2.11 

2.11 

2.11 

0.82 

0.82 

0.82 

0.82 

17792 

17792 

17792 

17792 

3.50 

0.27 

12717 


2.11 

0.82 

17792 


E  25627 

7001 

10200 

9420 

11821 

10923 

10152 

10509 

11889 

11715 

11653 

E  =  26556 

11755 

13917 

14459 

15516 

14724 

14454 

15127 

16632 

16575 

14933 

3.0,b  14.2,b  5.6,b  6.9,b  2.2lb  3.1 ,b  0.5,b  -0.5,b  -0.8lb  -1.0lb  -0.5lb  1.0lb 

CASE  2-1  BEAM  2 


1 9.5'“  14.2lb 


039  0-87  0.87  |  2. 

14.45  1.9  1.9 

51531  AAJCn  AA  1AA 


Ilf 


20469  20469  9738 


0.43 
=  13.35  4 
49530  o 


7716 

7716 

7716 

7716 

7716 

7716 

2.11 

2.11 

2.11 

2.11 

2.11 

2.11 

0.82 

0.82 

0.82 

0.82 

0.82 

0.82 

13447 

13447 

13447 

13447 

13447 

,»  0.024" 


1  ^ 

Q.  0.0  5"  0  01 


ooir  i1009"  0.008"  ooor  0006  0005  0004 


t  --  29364 

6089 

8237 

3097 

3371 

3005 

3626 

7032 

6248 

5979 

6356 

E  =  29545 

11293 

13407 

10829 

3414 

2563 

10402 

11231 

11647 

10926 

10379 

1.6,b  137lb  4.7lb  jojlb  2.2lb  5, lb  1.1,b  -0.5lb  -1.0lb  -1.3lb  -0.6lb  1.3lb 

CASE  3-2  BEAM  2 


V- ' 4 


V V: 


>«* 


CASE  3-2 


C-7 


m 


wi 


m 


1 

Bj 


$ 


i 


I  1  nntiR"  0006"  0.004"  0.003"  0.003"  0.002"  0.002"  0.002"  0.001"  000r 
0.014"  0.013"  0008 


a.  INPUT 


t  75394 

10336 

l  -  40611 

12882 

t  t 

39.5lb  30.9'b 

63.2lb 


5  gib  0  -1.2,b  -1.3lb  -1.0lb  -0.9,b  -0.6lb  0.7lb 

b.  OUTPUT 

CASE  4-3  BEAM  1 


6.7">  |4  jib  5.1lb  5.7lh  1.8lb  l.4,b  0.1,b  -0.4lb  -0.5lb 

CASE  4-3  BEAM  2 


-0.2lb  0.5lb 


0.0  27  00104- 

1  l  °“77  00058"  0.0041"  0.0031"  0.0025"  0.0021"  0.0018’  0.0015'  00013"  00012” 

S  ?  ?  »  .  •  '  •  k 


i  alb  -1.5lb  2.5lb  2.6lb  l.4lb  0-7'b  °'b  -0'3'b  -°-3'b  03lb  •°-1'b  °-2'b 


GRID  1 


0.0060"  0.0056”  0.0048"  0.0040"  0.0032"  0.0026"  0.0022"  0.0019"  0.0016"  0.0014"  0.0012”  0.001  V 


0.7lb  0.3lb  0.3lb  0.7lb  0.4lb  0.1 lb  0lb  -0.2lb  -0.2lb  0.1 lb  0.1 lb 

GRID  2 

CASE  4-3 


